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CHAPTER  I  -  INTRODUCTION 

NEPTUNE'  is  Nuclear  Electric  Propulsion  Technology  Using  NERVA.  Its  main 
components  consist  of  the  NERVA  nuclear  thermal  rocket,  acting  as  the  power  source  in 
a  closed  cycle  system,  a  magnetohydrodynamic  (MHD)  electricity  generator,  a 
magnetoplasmadynamic  (MPD)  thruster,  and  the  associated  devices  to  remove  waste  heat 
from  the  system  (i  e.  radiator).  It  is  intended  as  a  step  towards  a  high  power  nuclear 
electric  propulsion  technology. 

The  competition  between  conventional  chemical  propulsion  technologies  and  less 
established,  yet  potentially  more  rewarding,  nuclear  thermal  rocket  engines  began  in  the 
late  sixties  and  continues  to  the  present.  So  far,  chemical  systems  have  carried  the  burden 
successfully.  They  have  met  the  needs  of  the  missions  they've  been  designed  to  support. 
They  have  b'  .  oven  safe,  reliable,  high-thrust  propulsion  options,  and  this  has  always 
been  good  enough. 

NEPTUNE  represents  a  new  class  in  the  competition  for  high-thrust  propulsion.  It 
introduces  electric  propulsion  for  applications  in  ambitious  high-thrust  missions. 
Immediately  the  question  arises  -  if  nuclear  thermal  propulsion  was  not  necessary,  why  is 
nuclear  electric  propulsion  necessary? 

The  answer  is  simple.  Though  the  space  race  stagnates  at  times,  the  future  will 
certainly  bring  with  it  more  ambitious  and  more  frequent  departures  from  earth.  The  past 
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has  shown  that  chemical  propulsion  is  a  good  enough  way  to  meet  new  goals,  but 
eventually  there  will  be  a  calling  for  the  "best"  way. 

Is  nuclear  electric  propulsion  the  best  way  to  propel  a  spacecraft?  It  is  unlikely. 
There  is  a  lot  of  physics  to  be  resolved  before  words  like  "best"  will  become  appropriate. 
For  now,  however,  it  is  sufficient  to  say  nuclear  electric  propulsion  is  a  better  way  than 
either  chemical  or  nuclear  thermal  propulsion. 

This  presently  unsupported  statement  will  hopefully  gain  support  in  this  thesis,  but 
the  basics  are  clear.  The  breaking  of  nuclear  bonds  generate  more  energy  than  the 
breaking  of  chemical  bonds.  More  energy  is  available  from  less  material  if  nuclear  energy 
is  harnessed.  If  ever  there  was  a  forum  where  mass  must  be  minimized,  it  is  space. 

Nuclear  energy  is  made  for  space. 

The  comparison  of  nuclear  thermal  to  nuclear  electric  technologies  must  focus  on 
the  better  way  to  use  nuclear  energy.  It  is  natural  to  first  consider  nuclear  thermal 
rockets.  Nuclear  thermal  rockets  represent  the  use  of  nuclear  power  to  imitate  chemical 
rockets  —  the  basic  difference  being  the  origin  of  the  energy. 

But  this  should  not  be  the  only  difference.  Nuclear  energy  manifests  itself  in  a 
very  different  manner  than  the  energy  available  in  chemical  rockets  and  it  is  made  for  the 
closed  cycle.  That  is  its  nature.  The  nuclear  fuel  does  not  become  the  coolant  that  carries 
the  energy  out  of  the  system.  The  coolant  can  be  cycled  past  the  source,  as  it  is  in 
ground-based  plants.  It  is  different  from  a  chemical  source  in  that  capacity,  and  it  creates 
the  opportunity  to  use  the  energy  in  a  different  way  -  the  most  effective  way. 
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It  is  unnecessary  to  use  a  nuclear  sourced  system  to  mimic  chemical  systems.  The 
energy  should  be  used  in  the  best  way  to  apply  a  large  acceleration  to  a  reasonable  mass. 
The  following  chapters  will  explain  why  electromagnetic  forces  are  the  best  way  to  use 
this  energy. 

1.1  -  PROBLEM  STATEMENT: 

This  thesis  will  deal  with  the  problem  of  using  the  NEPTUNE  propulsion  concept 
as  the  primary  propulsion  source  for  a  manned  Martian  excursion.  It  will  define  a  round 
trip  Martian  trajectory  of  less  than  two  years  based  on  the  delta-v  parameter  of  that 
mission.  It  will  derive  a  budget  for  the  mass  of  a  nuclear  electric  propulsion  system  to 
undertake  that  mission,  and  it  will  analyze  NEPTUNE'S  ability  to  accommodate  that 
budget. 

In  the  event  that  NEPTUNE,  in  its  purest  form,  is  unable  to  meet  the  budgetary 
constraints  on  its  mass,  suggested  modifications  to  the  system  will  be  presented  and 
analyzed  which  will  allow  it  to  meet  them. 


1.2  -  NATURE  OF  THE  PROBLEM 

Although  creating  a  propulsion  system  specific  to  a  mission  seems  to  add  an 
element  of  complexity  to  the  analysis,  it  actually  adds  system  constraints  that  make 
analysis  possible.  This  strategy  allows  the  problem  to  be  approached  from  two  directions. 
First,  the  mission  is  known.  A  round-trip  mission  to  Mars  is  the  goal.  Moreover,  it  is  a 
manned  venture,  putting  definite  constraints  (less  than  two  years)  on  the  mission  length. 
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There  are  characteristics  of  a  round-trip  Martian  trajectory,  lasting  less  than  two 
years,  that  are  independent  of  the  propulsion  system  used.  Therefore,  NEPTUNE  must 
first  conform  to  these  constraints. 

The  problem  may  be  attacked  from  a  second  direction  because  the  propulsion 
system  to  be  used  in  the  mission  is  also  known.  NEPTUNE  has  definite  limitations.  The 
NERVA  source  is  limited  in  its  capacity  to  produce  power.  Its  electrical  generator  is 
limited  in  its  capacity  to  convert  power,  and  its  electric  thruster  is  limited  in  its  ability  to 
use  the  power.  Enough  is  known  about  NEPTUNE's  components  to  dictate  mass 
estimates  for  the  system  based  on  the  mission  demands.  However,  enough  is  left  unknown 
about  NEPTUNE's  components  to  allow  for  a  manipulation  to  accommodate  an  ambitious 
mission. 

NEPTUNE  is  a  nuclear  electric  propulsion  system.  Nuclear  electric  propulsion 
(NEP)  concepts  have  a  number  of  advantages  over  traditional  chemical  and  even  nuclear 
thermal  systems,  or  other  electric  propulsion  systems.  Because  they  use  electromagnetic 
field  forces  to  accelerate  propellants,  extremely  high  exhaust  velocities  are  possible. 
Additionally,  relatively  few  moving  parts  are  involved  in  thrust  production,  so  a  definite 
reliability  is  inherent.  Most  importantly,  however,  NEP  systems  bridge  the  gap  between 
electric  propulsion  and  thermal  systems  because  of  the  high  power  source  that  anchors 
them.  The  power  going  into  the  thruster  can  be  higher  than  otherwise  sourced  electric 
propulsion  systems  (i.e.  battery,  solar)  and  can  therefore  accelerate  more  propellant.  As  a 
result,  the  range  of  their  usefulness  might  not  exclude  higher  thrust  missions  if  sufficiently 
high  power  nuclear  sources  are  available. 
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The  manned  Mars  mission  from  low  earth  orbit  (LEO)  is  one  of  a  class  of  missions 
in  the  intermediate  thrust  range.  It  does  not  originate  from  the  earth's  surface,  so  it  does 
not  require  so  much  instantaneous  thrust  as  missions  having  to  overcome  surface  gravity. 
However,  the  "manned"  aspect  of  the  mission  limits  the  allowable  time  for  the  mission, 
which  leads  to  significant  thrust  requirements. 

This  thesis  will  analyze  NEPTUNE'S  ability  to  fulfill  the  requirements  of  this  class 
of  mission.  In  terms  of  specific  impulse,  it  operates  as  an  electric  propellant  accelerator, 
capable  of  exhaust  velocities  above  50,000  m/s.  In  terms  of  thrust,  it  has  a  source 
powerful  enough  to  produce  multimegawatts  of  power.  NEPTUNE  would  seem  to 
accommodate  that  intermediate  region  (which  encompasses  most  manned  planetary  space 
transfers)  in  which  the  craft  is  moved  only  from  planetary  orbit  to  planetary  orbit. 


CHAPTER  U  -  BACKGROUND 

Although  the  merits  of  NEPTUNE  are  fully  understood  only  when  the  system  is 
considered  as  a  whole,  it  is  useful  to  first  describe  the  elements  of  the  system  as  they  stand 
alone  to  provide  sufficient  background  as  to  how  they  succeed  together. 

2.1  -  THE  ELEMENTS  OF  NEPTUNE: 

REACTOR  SOURCE: 

The  subsystem  that  anchors  the  NEPTUNE  concept  is  a  compact,  high  power 
nuclear  rocket.  As  the  name  implies,  the  concept  would  ideally  employ  a  NERVA-style 
nuclear  thermal  rocket  as  the  power  source.  It  is  a  tested  reactor  and  rocket  engine 
combination,  projected  to  be  able  to  generate  flow  stagnation  temperatures  of  2500  K, 
thrust  levels  up  to  75000  lbs,  and  specific  impulse  values  above  750  seconds.2 

The  NERVA  reactor  is  of  prismatic  type.  In  it  a  gaseous  hydrogen  coolant  flows 
longitudinally  through  eight  canals  surrounded  by  a  homogeneous  mixture  of  uranium  fuel 
and  moderator.  As  the  hydrogen  gas  passes  over  the  fissioning  fuel,  the  heat  byproduct  of 
the  nuclear  reaction  taking  place  is  transferred  to  the  hydrogen  working  fluid. 

If  the  engine  is  acting  in  its  intended  capacity,  as  a  nuclear  thermal  rocket,  the  hot 
working  fluid  then  passes  through  a  nozzle,  transferring  the  heat  energy  into  kinetic 
energy  of  flow,  and  is  expelled  out  of  the  system.  The  momentum  transfer  of  the  high 
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velocity  hydrogen  produces  thrust  just  as  the  expulsion  of  the  byproducts  of  chemical 
propellants  do  when  they  are  thrown  from  traditional  chemical  systems.  Because  the 
molecular  mass  of  H2  is  low,  the  specific  impulse  of  this  type  of  system  is  high.  However, 
the  momentum  flux  (i.e.  thrust)  of  the  flowing  gas  is  reasonably  high  as  well  because  the 
reactor  is  able  to  add  heat  energy  to  a  substantial  total  mass  flow  of  hydrogen  gas. 


Figure  2.1  -  Prismatic  Reactor3 

Basically,  NEPTUNE  has  two  requirements  in  its  power  source.  It  must 
have  the  capability  to  transfer  heat  energy  at  the  rate  of  100  to  1000  MW  to  a  gaseous 
coolant,  and  secondly,  it  must  be  compact  enough  for  space  travel.  NERVA  was  designed 
with  these  very  same  goals  in  mind,  making  it  a  natural  selection  for  the  NEPTUNE 
system. 
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MHD  GENERATOR: 

Normally,  a  NERVA  nuclear  rocket  would  expel  its  high  energy  propellant 
through  a  nozzle  to  produce  thrust  in  an  open  thermodynamic  cycle.  However,  in 
NEPTUNE,  the  rocket  must  be  operated  on  a  closed-cycle.  The  gaseous  coolant  exits 
NERVA's  nozzle  and  enters  a  magnetohydrodynamic  (MHD)  generator. 

Magnetohydrodynamic  power  generation  is  the  conversion  of  the  kinetic  energy 
carried  by  a  moving,  electrically  conducting  fluid  (e.g.  plasma)  into  electrical  energy  by 
interaction  of  the  fluid  with  a  magnetic  field.  The  velocity  and  the  conductivity  each  tell 
something  about  the  gas  in  an  MHD  system.  The  velocity  indicates  how  much  energy  the 
gas  has.  The  conductivity,  c,  is  one  of  the  parameters  used  to  determine  the  magnetic 
Reynolds  Number: 

Rem=  opuL  (2.1) 

The  magnetic  Reynolds  Number  is  a  parameter  that  helps  to  define  the  ease  with  which  the 
energy  in  the  fluid  can  be  extracted  in  the  form  of  electric  power.  As  a  result,  the  most 
efficient  system  will  involve  a  gas  with  not  only  high  velocity,  but  high  conductivity.4 

The  drawing  of  the  fluid  element  uses  vector  geometry  to  describe  the 
electromagnetic  phenomenon  that  will  be  encountered  in  an  MHD  generator.  In  the 
volume  element,  consider  a  flowing  plasma  gas  with  velocity,  u,  and  scalar  conductivity, 
a.  As  such  a  gas  passes  through  a  magnetic  field  of  strength  B,  a  current  density  is 
induced  according  to  the  formula: 

/  =  c(£  +  u  x  (2.2) 

j :  current  density 
o :  conductivity 
u  :  flow  velocity 


B  :  magnetic  field 
E  :  electric  field 
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Figure  2.2  —  MHD  Fields5 

where  E  is  the  internal  electric  field  induced  to  support  the  load.  If  a  parameter  K  is 
defined  as: 

(2  3) 

then  the  current  density  can  be  rewritten: 

’  =  (1  -K)guB  (2.4) 

The  intent  of  an  MHD  generator  is  to  create  power  to  support  a  load.  The  amount  of 

power  it  can  produce  is  related  directly  to  the  current  density  via  the  formula: 

Po  =jE  =  K(\  -  K)ou2B2  (2.5) 

P0 :  generator  power  density 

If  the  above  relations  are  rearranged,  the  output  electric  field  and  power  per  unit  volume 

as  a  function  of  current  are  defined  as: 

E  =  uB-i  ^  (2.6) 

Po  =juB  -  (2.7) 

Therefore,  the  voltage  varies  linearly  and  the  power  quadratically  with  current  (as  is  the 

case  with  most  power  sources).4 


The  vector  diagram  pictures  a  jxB  force  acting  in  a  direction  opposite  the  gas  flow: 
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F=jB  =  (]  -K)auB2  (2.8) 

This  force  slows  the  gas  because  it  acts  opposite  to  the  direction  of  motion.  It  requires 

only  a  basic  understanding  of  the  law  of  conservation  of  energy  to  see  that  as  the  gas 

slows  and  the  corresponding  kinetic  energy  decreases,  the  energy  must  manifest  itself  in 

some  other  form,  rather  than  just  disappear.  In  the  case  of  the  MHD  generator,  it  is 

manifested  in  the  form  of  an  electric  field. 

Having  analyzed  the  electromagnetic  events  which  are  taking  place  in  an  MHD 
generator,  the  next  step  is  to  consider  the  device  as  a  working  system.  A  basic  MHD 
generator  system  can  be  broken  into  three  regions.  The  first  of  these  expands  a 
high-temperature  gas  through  a  nozzle  to  a  high  velocity.  The  high-temperature  gases 
may  be  produced  by  a  gaseous  fission  or  fusion  reactor  or  by  a  high-energy  combustor. 
From  there,  the  necessary  level  of  ionization  is  added  to  the  gas.  To  increase  the 
electrical  conductivity  of  the  gases,  they  are  often  seeded  with  an  easily  ionizable 
substance,  such  as  potassium  or  cesium.  This  region  may  or  may  not  be  included  with  the 
first.  It  will  depend  upon  the  amount  of  energy  initially  in  the  gas  and  the  difficulty  with 
which  the  gas  achieves  ionization. 

Finally,  the  gas  enters  into  a  third  region  where  the  actual  power  generation  takes 
place.  In  this  region  the  flowing  plasma  is  passed  through  a  magnetic  field  region  which 
contains  a  series  of  pairs  of  electrodes  at  either  side  of  the  magnetic  channel.  The 
electrode  pairs  face  each  other  and  lie  in  the  plane  parallel  to  the  plane  described  by  the 
flow  and  magnetic  field  directions,  which  are  at  right  angles  to  each  other. 

When  the  conducting  fluid  enters  the  magnetic  field,  interaction  will  occur  and  an 
induced  electric  field  is  established  normal  to  the  flow  and  magnetic  field  directions  and 
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between  the  electrode  pairs.  This  field  can  then  be  used  to  drive  a  current  via  the 
electrodes  to  an  external  load  where  it  is  allowed  to  perform  work.  In  doing  so,  the 
kinetic  or  potential  energy  of  a  gas  is  converted  into  electrical  energy. 

Consider  a  flowing  gaseous  conductor  in  a  uniform  magnetic  field  B.  An  e.m.f. 
(uxB)  is  set  up  at  right  angles  to  the  velocity  vector  and  the  field  as  shown  in  the  previous 
diagram.  If  the  Hall  number  (the  tendency  of  the  electrons  to  carry  current  in  the  direction 
of  the  flow)  is  very  small,  then  a  current  will  flow  parallel  to  the  induced  e.m.f.  The 
induced  e.m.f.  is  the  charged  particle  separating  force  caused  by  the  magnetic  field. 

r  \ 


emf  uxB 


Figure  2.3  —  Induced  E.M.F.4 

In  the  absence  of  electrodes,  as  in  the  above  diagram,  open  circuit  conditions  prevail  and 
the  electrons  will  flow  downwards  and  the  positive  particles  will  flow  upwards.  As  the 
charges  separate,  there  results  a  distribution  of  charge  (electrostatic  field)  which  opposes 
the  induced  e.m.f.  and  eventually  is  sufficient  to  balance  the  tendency  for  charge 
separation.  Therefore,  e.m.f.  and  further  current  flow  is  prevented,  E  =  -uxB. 
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In  order  to  allow  the  e.m.f.  to  drive  a  current  through  an  external  load,  it  is 
necessary  to  reduce  the  strength  of  the  electrostatic  field.  This  is  achieved  by  immersing 
electrodes,  which  are  connected  to  a  load,  into  the  fluid  as  shown  in  the  second  diagram. 
A  restricted  flow  of  electrons  then  occurs  from  the  cathode  through  the  fluid  to  the  anode 
returning  to  the  cathode  via  the  external  load.  The  driving  force  on  the  electrons  is  equal 
to  the  difference  between  the  induced  e.m.f.  and  the  electrostatic  field  £  and  is  called  the 
total  e.m.f.  The  actual  value  of  the  total  e.m.f.  depends  on  the  value  of  E,  which  depends, 
by  Ohm's  law,  on  the  total  resistance  of  the  electron  circuit.5 


Figure  2.4  -  Schematic  diagram  of  MHD4 
In  principle,  MHD  converters  draw  a  simple  analogy  to  conventional  electric 
generators.  A  rotating-type  generator  produces  electric  power  by  moving  a  metallic 
conductor  across  a  magnetic  field.  In  a  plasmadynamic  converter,  a  conducting  gas  moves 
through  a  magnetic  field.  An  electric  potential  difference  develops  across  two  electrodes 
which  are  placed  perpendicular  to  the  flow  direction  and  to  the  magnetic  field.  The  power 
output  of  the  generator,  like  that  of  the  conventional  rotating-type  generator,  is  a  function 
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of  the  magnetic  field  strength,  of  the  velocity  of  the  conducting  gas  through  the  field,  of 
the  linear  dimensions  of  the  device,  and  of  the  conductivity  of  the  working  medium  .4 

Ideally,  all  of  the  energy  taken  from  the  slowing  down  of  the  flow  would  be 
transferred  to  electrical  energy  manifested  in  MHD  power  to  drive  the  load.  In  reality, 
energy  must  go  to  generating  the  magnetic  field  and  powering  an  ionizing  source  to  create 
the  working  plasma.  This  adds  new  loads  to  the  circuit,  in  addition  to  the  thruster. 

Energy  is  also  lost  to  resistive  heating  stemming  from  obstructions  in  the  plasma  as  the 
current  moves  through. 

In  simple  terms,  an  MHD  generator  is  little  more  than  a  plasma  accelerator  acting 
in  reverse.  In  the  generator,  the  inlet  velocity  of  the  gas  is  greater  than  E/B.  In  this 
situation,  the  flow  will  be  decelerated,  and  current  will  be  driven  through  an  external  load 
by  the  motional  electromotive  force.  Reversing  the  conditions  would  accelerate  the  flow. 

Because  the  MHD  and  magnetoplasmadynamic  (MPD)  accelerators  are  so  similar, 
they  potentially  mate  into  attractive  and  complementary  systems  for  electric  space 
propulsion.  Since  the  spacecraft  must  already  be  set  up  to  handle  plasma  (possibly  seeded 
to  aid  ionization)  flow  through  its  channels  if  it  uses  an  MPD  accelerator,  little 
modification  will  be  necessary  to  accommodate  the  MHD  generator  to  support  the 
thruster.  This  kind  of  component  coupling  is  an  example  of  how  careful  engineering 
eliminates  complexity  and  the  associated  system  mass.  This  same  combination  of 
components  is  considered  in  the  NEPTUNE  concept. 
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MPD  THRUSTER: 


Figure  2.5  ~  MPD  thruster 


The  third  major  subsystem  associated  with  NEPTUNE  is  the  thruster.  A 
magnetoplasmadynamic  (MPD)  thruster  is  an  electric  thruster  which  employs  the  same 
cross-field  electromagnetic  theory  used  by  the  MHD  generator.  However,  the  thruster 
reverses  the  processes  and  takes  electric  energy  from  a  DC  source  and  uses  it  to  accelerate 
plasma  electromagnetically  for  the  purpose  of  producing  thrust. 


Figure  2.6  —  Ion/electron  reactions  to  Lorentz  force 
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MPD  thrusters  hold  a  major  advantage  over  other  types  of  electric  thrusters 

because  they  use  the  Lorentz  force  to  produce  thrust: 

F thrust  =jxB  (2.9) 

In  so  doing,  the  plasma,  with  its  component  electrons  and  ions,  though  responding  to 

fields  oppositely,  are  accelerated  in  the  same  direction  and  exit  the  craft  together.  Though 

the  gas  has  free  charges,  it  is  electrically  neutral  and  problems  associated  with  charge 

build-up  are  nonexistent. 

In  principle,  MPD  thruster  performance  should  improve  with  increased  power 
because  thrust  density  and  specific  impulse  are  increasing  functions  of  arc  current6: 

~F thrust  -  jf|(ln  77  +  J2  (2.10) 

J :  thruster  current 
r, :  anode  radius 
rc :  cathode  radius 
p0 :  permeability  of  free  space 


2.2  -  NEPTUNE  AS  A  SYSTEM: 

Chemical  propulsion  systems  have  traditionally  carried  the  burden  of  high  thrust 
space  ventures  because  they  are  able  to  expel  a  large  amount  of  mass  at  a  reasonable 
exhaust  velocity.  Consequently,  the  change  in  momentum,  the  thrust  force,  is  high  enough 
to  kick  spacecraft  in  the  right  direction  to  its  destination,  despite  the  forces  of  gravity 
resisting  that  change.  Electromagnetic  field  based  propulsion  systems  are  traditionally 
designed  for  less  ambitious  missions.  In  the  past,  the  power  contained  in  the  field 
accelerator  has  been  too  small  to  accelerate  enough  propellant  to  produce  high  thrust. 
What  makes  this  a  shame  is  the  fact  that  field  thrusters  can  accelerate  propellants  to  much 
higher  velocities  than  conventional  thermodynamic  accelerators  on  chemical  systems.  If 
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these  fields  were  only  powerful  enough  to  apply  their  force  to  a  larger  amount  of 
propellant,  the  momentum  changes  they  could  create  would  quickly  establish  them  as  the 
more  desirable  system,  in  terms  of  efficiency. 


Figure  2.7  -  NEPTUNE 


NEPTUNE  calls  for  the  application  of  the  theory  of  electric  acceleration  of 
propellants  to  a  high  thrust  impulsive  mission.  Just  as  the  ally  to  chemical  thrust 
production  is  high  mass  output,  the  ally  to  the  electrical  system  must  be  high  velocity 
output,  coupled  with  just  enough  mass  output  to  provide  the  required  amount  of  thrust. 
Even  at  its  best,  NEPTUNE  will  not  be  suited  for  truly  impulsive  missions.  The  thrust  is 
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not  applied  impulsively,  or  very  nearly  so,  as  is  the  case  with  conventional  chemical 
schemes.  Due  to  the  relatively  small  mass  output  and  despite  the  high  velocity  of  this 
output,  the  electric  thruster  will  only  be  capable  of  a  reduced  amount  of  thrust  and  if  it 
were  only  allowed  an  instant  of  firing,  its  effect  on  the  spacecraft  trajectory  would  be 
negligible.  But  how  reduced  must  this  thrust  level  be?  If  we  are  able  to  dump  more 
power  into  the  fields  of  acceleration  potential  than  they  generally  encounter,  can  we 
approach  the  momentum  change  of  a  chemical  system?  If  we  apply  the  thrust  at  a 
constant,  reasonably  high  level  for  a  longer  period  -  a  period  just  short  enough  to  make 
the  firing  time  look  impulsive  when  compared  to  the  extended  transfer  time  of  its  Mars 
mission  —  can  NEPTUNE  rival  chemical  systems?  This  section  describes  the  unique 
combination  of  NEPTUNE  components  that  will  allow  an  electric  propulsion  system  to 
operate  in  a  high  thrust  impulsive  domain. 

The  thrust  for  NEPTUNE  is  provided  by  a  multimegawatt  MPD  thruster,  capable 
of  generating  fields  with  enough  strength  to  accelerate  a  large  mass  flow  of  propellant. 
The  last  section  described  how  the  MPD  thruster  operates.  It  involves  no  moving  parts  in 
its  purest  sense  and  avoids  any  complications.  Since  complexity  often  leads  to 
inefficiency,  it  would  be  useful  to  keep  the  thruster  in  its  purest  state  in  NEPTUNE. 

Since,  the  electric  field  of  an  MPD  thruster  is  a  static  field  simply  running  from  an  annular 
anode  to  a  cylindrical  cathode  it  needs  simple  DC  power  to  create  a  difference  in  potential 
across  its  curved  "plates".  However,  to  be  useful  for  the  NEPTUNE  mission,  the  fields  in 
the  thruster  must  be  powerful  and  the  source  to  drive  them  must  be  likewise.  Therefore, 
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NEPTUNE  needs  a  device  to  supply  a  large  magnitude  of  DC  power  to  mate  with  its 
unique  MPD  thruster. 

This  leads  us  a  step  backward  along  NEPTUNE  to  the  MHD  generator.  The 
MHD  generator  electromagnetically  decelerates  a  high  velocity  conducting  fluid  and 
converts  the  flow  energy  to  direct  current  electrical  energy.  Obviously,  if  the  amount  of 
energy  contained  in  the  flow  is  sufficient,  this  DC  source  will  lead  perfectly  into  the  MPD 
thruster.  Absent  are  the  converters,  the  moving  parts,  and  the  inefficiencies. 

The  key,  however,  is  to  ensure  that  enough  DC  power  is  coming  from  the  MHD 
generator.  So  we  take  yet  another  step  backward  along  the  NEPTUNE  cycle.  What  is 
needed  is  a  device  that  efficiently  creates  a  very  high  energy  flow  conducting  gas.  The 
conductivity  part  of  this  requirement  is  handled  by  an  arc  discharge  just  prior  to  the 
entrance  of  the  generator.  This  ionizes  the  previously  neutral  gas.  Additionally,  the  arc  is 
composed  of  a  DC  field  which  can  be  supported  by  the  same  DC  generator  which  it  helps 
to  run.  The  key  being  the  coupling  of  devices  that  work  naturally  together. 

The  second  part  of  the  requirement  is  for  a  high  velocity  gas.  Chemical  reactions 
can  result  in  high  temperature  gas  to  be  expanded  to  a  high  velocity,  but  the  magic  of  the 
chemical  reaction  is  not  something  that  can  be  repeated  in  the  closed  cycle.  It  calls  for 
something  more  simple  —  something  to  act  as  a  very  powerful  stovetop  burner,  so  as  the 
fluid  passes  by  it,  multimegawatt  heat  transfer  occurs.  That  is  what  a  nuclear  reactor 
does.  It  adds  energy  to  a  fluid  in  the  form  of  heat  and,  because  of  the  potential  of  the 
nuclear  bond,  the  magnitude  of  this  energy  is  very  high.  Furthermore,  by  using  the 
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NERVA  engine,  the  compact  reactor  is  included  to  supply  the  energy  in  the  form  of  heat 
and  a  nozzle  is  attached  to  transfer  the  heat  into  kinetic  energy  in  the  flow. 

The  NERVA  rocket's  intended  application  was  as  a  nuclear  thermal  rocket.  In  this 
class  of  propulsion  device,  propellant/coolant  is  taken  from  a  storage  tank,  passed  over 
the  very  hot  nuclear  reactor  source,  taking  with  it  the  source's  heat,  and  finally  accelerated 
with  a  nozzle  which  converts  the  heat  energy  in  the  gas  to  kinetic  energy  of  the  flow. 

Now  the  flow  of  this  high  velocity  fluid,  normally  expelled  to  the  environment  and 
leading  directly  to  thrust,  may  be  contained  and  used  on  the  closed  cycle.  It  turns  out  that 
this  is  just  the  type  of  very  high  energy  content  gas  that  can  be  exploited  by  an  MHD 
generator  to  create  very  high  levels  of  electrical  energy.  This  can,  in  turn,  be  exploited  by 
the  MPD  thruster,  using  large  electric  fields  to  expel  moderate  amounts  of  propellant  at 
very  large  velocities  and  producing  thrust. 

NEPTUNE  as  a  propulsion  system  represents  the  combination  of  components 
described  above.  It  is  founded  on  three  basic  advantages  that  will  be  noted  throughout  the 
thesis.  The  first  is  the  incorporation  of  the  same  high  power  producing  technologies  that 
are  performing  on  the  ground,  condensed  for  use  in  space.  The  second  and,  in  my 
opinion,  most  significant  is  the  manipulation  of  field  theory  for  our  own  benefit.  This 
represents  the  manipulation  of  nature,  or  the  use  of  force  in  the  form  nature  creates  it. 
Finally,  NEPTUNE  is  about  efficiency  --  simply  assembling  pieces  that  belong  together 
and  creating  success  through  simplicity. 
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However,  NEPTUNE  is  only  one  combination  of  these  ideas.  They  are  by  no 
means  limited  to  NEPTUNE.  If  not  in  NEPTUNE,  these  ideas  will  surface  in  whatever 
propulsion  system  develops  to  support  the  needs  of  the  future. 

NEPTUNE'S  us*  of  a  nuclear  source  represents  how  h  exploits  the  first  advantage. 
Whether  or  not  NERVA  is  the  source  of  choice,  the  harnessing  of  nuclear  energy  is  a  clear 
improvement  over  the  use  of  chemical  energy.  The  forces  holding  the  nucleus  of  an  atom 
in  place  are  simply  stronger  than  the  bonds  between  chemicals.  As  a  result,  it  makes  sense 
that  more  energy  can  be  derived  from  the  process  of  splitting  the  nucleus,  than  the 
molecule.  Energy  density  is  key  in  space  because  mass  is  such  an  overriding  concern.  The 
same  bulky  systems  which  have  great  success  on  land  do  not  translate  to  the  realm  of 
space.  While  high  power  is  still  important,  low  mass  is  key.  Nuclear  power  represents 
very  concentrated  energy.  For  this  reason,  it  belongs  in  space. 

The  second  advantage  is  more  abstract.  It  involves  a  subject  that  is  easily  the  most 
poorly  understood  by  the  scientific  community.  Field  forces  represent  action  at  a  distance. 
They  are  excited  by  mass,  in  the  case  of  gravity,  or  charge,  in  the  case  of 
electromagnetism.  The  final  understanding  of  the  relationship  between  mass  and  charge 
will  be  the  final  understanding  of  matter  and  the  universe,  but  what  is  important  to  know 
in  the  interim  is  that  nature  has  organized  the  universe  with  the  use  of  fields.  They  are  a 
natural  occurrence.  They  are  how  nature  wants  energy  to  exist. 

In  NEPTUNE,  the  NERVA  rocket  is  used  very  differently  than  intended,  to  create 
a  field  based  system.  A  closed  cycle  of  gas  is  the  platform  in  which  energy  is  carried.  It  is 
loaded  on  at  the  reactor  in  the  form  of  heat.  It  is  rearranged  by  the  nozzle  into  kinetic 
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energy  of  flow.  Then  it  is  off-loaded  in  the  MHD  generator  as  electricity,  only  to  return  to 
the  reactor  source  and  complete  the  cycle.  This  energy  then  leads  to  another  device  that 
uses  electromagnetic  field  forces  to  accelerate  an  ionized  gas  (plasma)  from  the  system  to 
produce  thrust.  Because  nature  must  undergo  no  conversion  when  fields  are  used, 
electromagnetic  forces  can  attain  much  higher  exhaust  velocities  than  are  possible  by 
expanding  a  fluid  thermodynamically,  as  in  nuclear  thermal  or  chemical  systems. 

The  magic  of  field  forces  is  key  to  systems  of  the  future  because  the  understanding 
that  will  come  with  time  seems  to  lead  to  unheard  of  potential.  As  stated  earlier,  the 
breakthrough  that  explains  why  they  can  act  as  they  do  is  the  breakthrough  that  explains 
what  makes  up  our  surroundings. 

The  final  key  to  NEPTUNE  is  basic,  but  important  because  of  its  range  of 
usefulness.  This  is  the  idea  that  waste  may  be  removed  from  a  system  if  the  time  is  taken 
to  create  a  system  from  components  that  belong  together.  Absent  from  NEPTUNE  are 
the  moving  parts  characteristic  of  most  propulsion  systems.  This  translates  into  system 
efficiency  and  reliability.  But  even  more  importantly,  what  is  absent  from  NEPTUNE  are 
unnecessary  energy  conversion  devices  and  the  reductions  in  efficiency  they  cany  with 
them.  Because  a  high  thrust  MPD  device  requires  a  powerful  DC  source,  an  MHD 
generator  is  used.  Because  the  generator  needs  high  speed  flow,  the  NERVA  engine  is 
used.  The  moral  of  the  story  is  that  maximum  efficiency  is  achieved  through  optimum 
component  matching. 


CHAPTER  HI  -  APPROACH 


The  intent  of  this  chapter  is  to  present,  in  detail,  the  methodology  used  in 
determining  NEPTUNE's  ability  to  support  a  manned  Mars  venture.  Although  this 
chapter  will  not  contain  any  of  the  actual  calculations  or  results,  the  exact  processes  to 
obtain  the  results  will  be  explained  in  detail. 

The  methods  used  to  evaluate  NEPTUNE's  ability  to  support  a  manned  Mars 
mission  can  be  broken  into  three  parts,  each  explained  in  a  separate  section  of  this  chapter. 
The  ultimate  goal  is  to  define  a  maximum  mass  for  the  NEPTUNE  propulsion  concept  if  it 
is  to  be  expected  to  carry  out  a  given  round  trip  Mars  mission  scenario.  In  order  to  define 
this  upper  limit,  an  equation  to  relate  the  propulsion  system  mass  to  parameters  describing 
the  mission  is  necessary.  This  equation  is  the  subject  of  the  first  section  of  this  chapter. 

The  next  step  is  to  calculate  the  key  parameter,  used  in  the  mass  equation, 
describing  the  mission,  based  on  the  allowable  mission  time  of  flight.  This  parameter  is  the 
delta-v  and  it  represents  the  magnitude  of  the  velocity  change  required  to  carry  out  the 
Mars  transfer.  The  method  by  which  this  parameter  is  determined  will  be  the  subject  of 
the  second  part  of  this  chapter. 

Finally,  the  third  section  will  derive  an  equation  to  analyze  the  elements  of 
NEPTUNE  and  evaluate  their  ability  to  meet  the  mass  budget  defined  by  combining  the 
first  two  sections  of  the  chapter.  In  other  words,  it  will  compare  what  must  be  possible 
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for  the  mission  to  succeed  to  what  is  actually  possible  based  on  the  components  which  are 
NEPTUNE. 

The  methods  are  described  in  a  sensible  order  in  this  chapter  —  starting  with  the 
primary  unknown  for  NEPTUNE,  the  system  mass,  and  moving  to  those  parameters 
necessary  to  solve  for  this  unknown.  Howe’  er,  in  the  actual  calculation  process,  it  is 
more  convenient,  in  order  to  test  every  possible  scenario,  to  use  a  different  order.  More 
specifically,  the  following  process  outline  will  guide  the  calculations  for  the  first  part  of 
this  thesis: 

1  —  calculate  the  delta-v  required  for  the  mission 

2  —  solve  for  the  specific  mass  of  a  system  with  the  capabilities  of  NEPTUNE 

to  perform  the  mission 

3  —  compare  this  specific  mass  value  to  another  compiled  from  the  specific 

masses  of  the  component  parts  of  NEPTUNE 

In  performing  this  process,  we  may  compare  how  "good"  NEPTUNE  must  be  to  how 
"good"  NEPTUNE  is.  The  results  from  this  analysis  will  answer  the  question  presented  in 
the  problem  statement  —  is  NEPTUNE  able  to  support  a  manned  Mars  venture? 


3.1  -  SPECIFIC  MASS  BUDGET: 

The  previous  section  draws  a  relationship  between  mission  time  of  flight  to  delta-v 
for  the  mission.  Now  a  relationship  to  propulsion  system  mass  must  be  added.  This  basis 
for  the  relation  of  delta-v,  and  therefore,  time  of  flight,  to  propulsion  system  mass  is  the 
rocket  equation7: 

1 Jo  _ 


n^  :  initial  mass  of  spacecraft 
mf :  spacecraft  mass  after  firing 


(31) 
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Av  :  change  in  velocity  required 
ue  :  exhaust  velocity  of  propellant 

However,  in  NEP  systems,  this  equation  holds  different  meaning  than  in  traditional 
chemical  systems.  In  chemical  engines,  the  difference  between  the  initial  mass  and  the 
final  mass  is  the  mass  of  the  propellant  and  propulsion  system,  since  the  propellant 
dominates  the  mass  of  the  propulsion  system.  In  NEPTUNE,  however,  the  propulsion 
system  mass  cannot  be  neglected  compared  to  the  mass  of  the  propellant.  In  other  words, 
the  initial  mass  can  be  defined: 

^propdimt  ®propuIiion  iy«tem  «V»d  (32) 

and  the  final  mass. 

mf  ^propulsion  «yjtero  +  «Uw  (33) 

This  leads  to  a  new  rocket  equation,  designed  for  concepts  with  power  sources  separate 
from  propellant7: 

<3-4> 

Consider  a  vehicle  that  will  travel  to  Mars  in  a  short  time.  If  the  amount  of  thrust 
it  is  capable  to  apply  is  defined,  it  would  like  to  maintain  the  vehicle  weight  at  the  lowest 
value  possible  so  that  thrust  force  will  have  more  of  an  effect  in  accelerating  it.  In  other 
words,  it  is  important  to  attain  a  high  thrust-to-weight  ratio.  The  thrust  will  ideally  be 
high  and  the  weight  low.  It  is  obviously  an  advantage  if  the  amount  of  propellant  used  to 
create  the  thrust  is  low.  In  order  for  that  to  happen,  the  velocity  at  which  the  propellant  is 
expelled  must  be  very  high  to  compensate.  This  describes  the  theory  behind  electric 
propulsion,  but  the  theory  is  complicated  by  other  considerations.  In  order  to  support  the 
electric  propulsion,  some  kind  of  electricity  generator  must  be  a  part  of  the  propulsion 
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system.  Obviously,  if  the  mass  of  this  electricity  production  plant  is  large,  it  will  negate 
the  savings  achieved  by  the  lower  propellant  mass  requirement. 

The  factor  of  interest  in  this  thesis  is  the  mass  of  the  propulsion  system.  It  is  what 
will  make  NEPTUNE  different  than  other  systems  designed  for  interplanetary  travel. 
Typically,  the  mass  of  the  electricity  plant  has  significantly  diminished  the  benefit  of  the 
electric  thruster.  However,  also  typically,  electric  propulsion  has  not  been  applied  to  high 
thrust  endeavors.  Fortunately,  scaling  to  these  high  thrust  and  high  power  missions  is  not 
linear  and  if  electric  thrusters  can  be  designed  to  handle  the  higher  power,  a  space 
electricity  generator  is  conceivable  with  low  enough  weight  to  make  the  resulting  vehicle  a 
definite  advantage  * 

Because  scaling  will  be  such  a  factor  in  the  analysis,  the  mass  of  the  propulsion 
system  is  broken  into  mass  per  unit  power,  or  the  specific  mass: 

^propuliion  lystm  ^prop  jyttP  prop  lyit  0-5) 

Ppopiyst :  specific  mass  of  propulsion  system 
Pprop.**  :  power  of  propulsion  system 

In  this  manner,  the  power  requirements  for  the  mission  can  simply  be  coupled  with  the 

specific  mass  to  lead  to  a  overall  propulsion  system  mass. 

If  the  substitution  is  made  for  the  propulsion  system  mass  in  the  external  source 

rocket  equation,  terms  can  be  rearranged  to  obtain: 

ntprop  =  (ntpoy  +  MVX*AW,“*  -  1)  (3.6) 

Over  the  life  of  the  propulsion  system,  it  will  be  expected  to  apply  a  total  impulse, 

or  change  in  velocity,  to  complete  the  mission.  This  will  equal  the  thrust  that  will  be 

produced,  throughout  the  entire  mission,  by  the  propulsion  system: 


I, :  total  impulse 


(3.7) 
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t  :  thruster  firing  time 
F :  thruster  force 

This  total  impulse  can  also  be  related  to  specific  impulse,  a  parameter  of  the  thruster  to  be 
used3: 

4  =  /vgoJo  md  (3.8) 

Iq, :  specific  impulse 

Note  that  the  mass  flow  rate  integrated  over  the  entire  time  of  burning  will  simply  be  the 
mass  of  the  propellant.  The  propellant  mass  is  seen  also  in  the  rocket  equation,  so  a  new 
equality  is  set  up: 

l,=Ispgo[(e*Wu-  -  1  Xmpay  +  MV)]  (3.9) 

Next,  if  we  consider  the  propulsion  system  in  its  basic  form  as  a  machine  that  is 

not  throttleable,  meaning  a  machine  that  applies  only  one  constant  level  of  thrust  whenever 

it  is  burning,  the  relationship  with  the  total  impulse  is  exploited  to  get  an  equation  for 

thrust: 

F,  =  ^[(f““--lXmw  +  P^„)]  (3  10) 

Finally,  since  the  power  of  the  thruster  is  related  to  its  thrust  force  by: 

P,  =  \FtUex  (3.11) 

a  substitution  into  the  thrust  equation  derived  above  returns  a  relation  for  the  thruster 

power: 

Pi  =  #((«Av/““  -  1)(»W  +  P  psPps)]  (3.12) 

This  power  relationship  is  very  useful  since  a  great  deal  is  known  about  the  power  source 

in  the  NEPTUNE  system.  NERVA  is  a  compact  nuclear  reactor  with  definite  capability 

limitations.  Its  maximum  power  level  is  established.  Therefore,  knowing  how  much 

power  is  produced  by  the  source,  a  fairly  strong  estimate  may  be  made  as  to  how  much 

power  will  eventually  reach  the  thruster  and  be  available  to  produce  thrust.  This  estimate 
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will  come  from  general  understandings  of  the  electricity  generator  and  the  electric  thruster 
and  their  associated  efficiencies: 

(3  13) 

Pmpd  =  empdeBnP>cUr«  (314) 

Having  derived  the  basic  equations  used,  the  procedure  to  calculate  a  specific  mass 
budget,  or  the  highest  total  specific  mass  allowable  for  the  NEPTUNE  propulsion  system 
to  complete  the  intended  mission,  combining  the  parameters  of  the  mission  with  the 
characteristics  of  the  propulsion  system  is  outlined  below: 


1  —  Define  the  available  power.  This  is  dependent  upon  the  power  production 

capabilities  of  the  NERVA  source.  It  is  a  known  parameter. 

2  -  Use  the  efficiencies  of  the  thruster  and  generator  subsystems  of  the  propulsion 

system  to  calculate  the  power  that  will  be  available  for  the  thruster. 

3  —  From  a  knowledge  of  the  exhaust  velocities  obtainable  by  the  thruster  to  be 

used,  calculate  the  required  mass  flow  rate  of  the  propellant. 

Pmpd  =  {  rrimpd  u2ex  (3.15) 

4  --  Determine  the  total  mass  of  the  propellant  by  combining  the  mass  flow  rate 

with  the  total  time  during  which  the  propulsion  system  will  be  applying  thrust. 
The  firing  time  is  defined  by  a  percentage  of  the  total  transfer  time,  to  maintain 
an  impulsive  trajectory.  It  does  not  take  into  account  capabilities  of 
NEPTUNE  to  operate  on  a  continuous  basis,  although  electric  thrust  systems 
are  generally  used  in  this  capacity. 

ffip—tfltbun  (3-16) 

5  —  Plug  values  into  the  external  source  rocket  equation  and  solve  for  the  specific 

mass  budget,  0,,. 

ZmpdSgtnP s  —  —  lytTlpay  -+  fiplP pi)] 


(317) 
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3.2  -  DELTA-V  FOR  IMPULSIVE  MARS  MISSION: 

In  any  impulsive  transfer  calculation  it  is  always  important  to  minimize  the  delta-v 
required  for  the  mission.  In  manned  transfers,  however,  it  is  just  as  important  to  minimize 
the  transfer  time  for  the  mission.  The  two  parameters  are  directly  related  and  clearly 
competing.  A  short  transfer  time  calls  for  a  more  direct  trajectory,  leaving  less  work  to 
the  forces  of  gravitational  attraction  between  the  bodies  and  more  for  the  propulsion 
system.  This  translates  directly  into  a  high  delta-v  requirement  for  the  mission  and  a  larger 
propulsion  system  or  longer  firing  time.  Since  the  firing  time  is  limited  by  the  assumption 
that  the  transfer  must  be  legitimately  considered  impulsive,  the  obtainable  delta-v  is 
limited,  in  large  part,  by  the  thrust  capability  of  the  propulsion  system. 

Different  mission  scenarios  are  examined  in  the  program  MARS  TRIP  (Appendix 
A).  The  mission  delta-v's  calculated  vary  with  the  transfer  time  and  planetary  geometry. 
Because  the  planetary  geometries  repeat  periodically  (particularly  when  considered  in 
two-dimensional  space),  the  program  is  set  up  to  evaluate  the  mission  at  the  entire  cycle  of 
geometries  for  a  given  transfer  time.  Transfer  time  limits  are  imposed  to  make  the  mission 
realistic  for  manned  travel  and  realistic  for  propulsion  system  capabilities.  The 
mathematical  technique  used  in  the  program  is  outlined  below. 

The  trajectory  analyzed  in  defining  the  Mars  transfer  is  basic.  The  mission  begins 
from  a  low  earth  parking  orbit.  From  this  orbit,  an  impulsive  bum  is  applied  to  set  the 
spacecraft  on  its  way  towards  Mars.  The  magnitude  of  the  thrust  is  dependent  on  the 
amount  of  time  allowed  for  this  outbound  transfer.  As  the  craft  approaches  the  Martian 
sphere  of  influence  on  a  trajectory  that  is  hyperbolic  with  respect  to  the  planet,  a  second 
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bum  is  used  to  insert  the  vehicle  into  a  closed  orbit  of  Mars.  During  the  time  that  the 
landing  party  is  on  the  planet's  surface,  the  transfer  craft  continues  to  orbit.  After  the  stay 
time  has  passed,  a  final  bum  sends  the  spacecraft  on  its  return  transfer. 

Obviously,  the  geometry  of  the  Sun-earth-Mars  system  plays  a  major  factor  in  the 
magnitude  of  the  delta-v  necessary  to  complete  the  round-trip  venture. 

Specifically  the  analysis  process  can  be  broken  down  as  follows: 

Delta-v  at  point  A: 

Point  A  is  defined  as  the  point  in  low  earth  orbit  (LEO)  of  initial  impulsive  bum. 
The  bum  is  performed  in  a  400  km  parking  orbit  of  earth.  The  desired  velocity  vector  is 
calculated  using  a  simple  relation: 

vM^m  =  5U/eeo  +  iW«-  (318> 

v,^  :  velocity  of  satellite  with  respect  to  the  sun 
IWo :  velocity  of  satellite  with  respect  to  earth 
v  .lim  :  velocity  of  earth  with  respect  to  the  sun 

All  velocities  may  be  considered  in  this  fashion  if  standard  units  are  used. 

In  the  above  relation,  vltL,vm  may  be  considered  the  velocity  of  the  spacecraft  at  the 
perihelion  of  its  transfer  orbit  to  Mars  in  an  heliocentric  or  sun-centered  reference  frame. 
This  point  will  always  be  the  perihelion  and  as  time  of  flight  restraints  dictate  more  direct 
missions,  the  size  of  the  elliptical  outbound  transfer  will  grow  so  that  Martian  orbit  radius 
is  reached  earlier  in  the  orbit.  This  vector  may  be  calculated  according  to  the  formula: 

sm(v,)P+^[«  +  cos(v,))£  (3.19) 

:  gravitational  parameter  for  the  sun 
p :  semi-latus  rectum 
vt :  true  anomaly  at  point  a 
e :  eccentricity 
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In  order  to  define  this  vector,  the  parameters  p  and  e  of  the  transfer  orbit  must  be  known. 
This  proves  to  be  a  profound  problem,  solved  using  an  iterative  process  in  subroutine 
ECCENT  (Appendix  A). 


Subroutine  ECCENT  addresses  the  problem  of  determining  the  size  of  an  orbit 
necessary  to  meet  limits  imposed  by  outbound  trip  times  to  Mars.  The  first  step  is  to 
assume  a  value  for  eccentricity  of  the  transfer  orbit,  e.  We  could  begin  with  a  value  of 
zero  for  eccentricity,  corresponding  to  a  circular  orbit.  However,  this  would  result  in 
more  iterations  than  necessary  and  a  longer  computer  run  time.  If  the  initial  eccentricity 
guess  is  instead  based  upon  a  Hohmann  transfer  eccentricity,  calculations  may  be  reduced. 
A  Hohmann  transfer  represents  the  lowest  energy  impulsive  trajectory  that  will  allow  the 
spacecraft  to  go  between  earth  and  Mars.  If  this  trajectory  were  followed,  the  earth 
would  be  at  the  point  of  perihelion  of  the  orbit,  while  Mars  would  be  at  the  apohelion 
point  just  as  the  spacecraft  arrives,  as  shown  in  the  diagram  of  a  Hohmann  transfer. 

This  orbit  has  a  definite  eccentricity: 


c  = 

(^t) 


e :  eccentricity 


(3.20) 
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r, :  distance  from  sun  to  point  a 
rb :  distance  from  sun  to  point  b 

Using  this  eccentricity  as  an  initial  guess,  the  semi-major  axis  is  found: 

r0 

a=M 

Also,  the  semi-latus  rectum  is  found: 
p  =  a(  1  -e2) 


(321) 


(3.22) 


Now,  find  the  angle  between  r,  and  the  radius  vector  for  the  spacecraft  when  it  reaches 
Mars  at  point  B  (in  this  case,  at  the  apohelion): 

cosv=mS>  <323) 

From  the  true  anomaly,  it  is  helpful  to  model  the  transfer  orbit  as  a  circular  orbit  and 

calculate  a  new  angle,  the  eccentric  anomaly,  based  on  this  model9: 

cos  (3.24) 

E :  eccentric  anomaly 

Now  a  time  of  flight  corresponding  to  the  eccentricity  guess  is  calculated: 


TOF  =  ;£[£  -  e  cos£]  (3.25) 

Finally,  this  new  time  of  flight  is  compared  to  the  allowed  time  of  flight.  If  they 
are  the  same,  this  is  the  orbit  to  use  for  the  transfer.  If  they  differ,  the  eccentricity  must  be 
increased  by  some  amount  to  make  the  transfer  orbit  larger  so  that  the  spacecraft  reaches 
Martian  radius  earlier  in  the  orbit  and  the  time  of  flight  between  the  planets  is  lower.  This 
subroutine  will  accept  orbit  parameters  as  correct  once  they  predict  the  time  of  flight 


within  half  of  one  earth  day. 


The  next  step  is  to  calculate  the  energy  necessary  to  create  a  transfer  orbit  of 
appropriate  size.  Obviously,  lower  trip  times  will  require  more  eccentric  parabolic  or 
hyperbolic  orbits  and  more  energy.  The  problem  is  broken  up  so  that  the  velocity  change 
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is  calculated  at  two  points.  The  first  is  the  point  where  the  transfer  orbit  begins  (the 
periapsis  of  the  orbit)  and  the  other  is  the  point  where  the  transfer  ellipse  interseas  the 
circular  orbit  of  Mars  Since  a  parabolic  orbit  retains  constant  energy  along  its  path, 
calculations  at  other  points  are  unnecessary. 

Now  the  parameters  for  the  outbound  Martian  transfer  are  known  and  the  vector 
velocity  that  the  spacecraft  must  have  at  earth  departure  in  order  to  undertake  this  transfer 
can  be  calculated.  The  coordinate  system  being  used  is  defined  by  the  plane  of  the  transfer 
orbit  from  earth  to  Mars.  The  P  direaion  is  defined  as  the  direaion  of  periapsis  in  this 
orbit,  which  is  also  a  position  veaor  to  pt.  A.  Consequently,  o,  =  0  radians. 

The  second  veaor  to  be  found  is  vftca.IulI.  The  same  formula  is  employed: 

=  JW  sin(v.)/>+  JW  co  5(v„)0  (3.26) 

In  the  previous  equation,  since  the  orbit  of  the  earth  about  the  sun  is  considered  circular,  p 
=  r,  and  e  =  0 

Now  these  veaors  are  amply  subtraaed  to  determine  the  velocity  that  the  satellite 
must  have  with  respea  to  earth  at  the  periapsis  point  of  its  heliocentric  transfer  orbit  to 
Mars. 


Originally,  the  satellite  will  have  a  velocity  at  point  A  of: 


rB :  radius  of  circular  orbit  about  earth 
va :  velocity  of  satellite  in  circular  earth  orbit 
Po :  gravitational  parameter  of  earth 

To  escape  the  influence  of  the  earth,  it  will  need  to  first  attain  escape  velocity,  but  this  will 
not  be  enough.  Escape  velocity  is  defined  as  the  exaa  velocity  necessary  for  the  craft  to 
escape  earth’s  gravity.  As  the  distance  from  the  earth  becomes  infinite,  the  velocity  of  the 
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craft  will  be  exactly  zero.  However,  for  our  purposes,  the  craft  must  have  a  velocity  of 
as  the  gravity  of  earth  is  overcome.  Consequently,  is  our  hyperbolic  excess 
velocity,  or  the  velocity  of  the  craft  when  the  distance  from  the  earth  to  the  craft  may  be 
considered  infinite.  We  may  calculate  that  velocity  necessary  in  LEO  to  give  us  y^^ 
outside  the  earth's  sphere  of  influence  by  the  following  relation. 

vL  =  vL  +  vL/*eo  (3-28) 

vto  :  required  bum-out  velocity 
v^ :  earth  escape  velocity  at  LEO 

V  sat /sun  ~  V sat/geo  "h  V  geo/sun  (3-29) 


Then  the  delta-v  that  will  be  necessary  is  a  simple  matter: 


&-a  ~  \-bo~~cs 


(3.30) 


Delta-v  at  point  B  (#1): 

For  the  Hohmann  transfer  problem,  this  delta-v  at  point  B  is  easily  calculated  much 
as  the  delta-v  at  A  was,  but  the  problem  is  more  complex  for  the  case  of  more  direct 
transfer  orbits.  This  is  because  the  vector  direction  of  the  transfer  orbit  velocity  at  point  B 
is  not  in  the  same  direction  as  the  circular  satellite,  or  Martian,  orbit  velocity  direction  for 
this  case.  Instead,  the  two  velocity  vectors  will  be  off  by  an  angle  6,  determined  by  the 
position  in  the  elliptical  transfer  orbit  where  the  spacecraft  encounters  Mars.  This  angle 
may  be  calculated  with  the  formula: 
cos6  =  ^vM 


(3.31) 
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Now,  the  thrust  must  not  only  be  applied  to  slow  down  and  circularize  the  orbit,  but  to 
manipulate  the  velocity  vector  into  the  proper  direction  for  a  circular  satellite. 

This  constitutes  the  second  bum  of  the  NEPTUNE  mission,  when  the  spacecraft 
encounters  the  vicinity  of  Mars.  In  order  to  determine  the  magnitude  of  this  delta-v,  it  is 
important  to  first  determine  the  exact  location  in  the  transfer  orbit  when  the  spacecraft 
reaches  Mars.  In  other  words,  find  the  true  anomaly  of  Mars  encounter.  The  true 
anomaly  is  the  angle  between  the  perihelion  direction  (in  this  case,  point  A)  and  the 
position  vector  to  the  spacecraft  at  any  point.  The  true  anomaly  at  point  B  can  be  found: 

^  =  (332) 

where  rb  is  the  magnitude  distance  from  the  sun  to  point  B,  or  Mars.  Now  the  position 
and  velocity  vectors  of  the  spacecraft,  at  the  point  of  Mars  encounter,  may  be  calculated 
using  the  familiar  relations: 

rb  =  r&cos  (\b)P  +  r  j,sin  (vb)Q  (3.33) 

£<*.</ sun)b  =  sin (yb)P  +  j?[e  +  cos(vj,)](?  (3 .34) 

Additionally,  the  velocity  vector,  v^^,  may  be  found  with  the  same  formula  as  above, 
assuming  p  is  replaced  with  rb  and  e  is  assumed  to  be  zero  since  the  Martian  orbit  is  nearly 
circular  about  the  sun. 

The  unknown  is  the  velocity  of  the  spacecraft  with  respect  to  Mars  at  point  B. 

This  may  be  found  with  simple  vector  addition: 


¥(M/MB)b  *  ^(«Mn]b  +  £(Mnfeii)b 


(3.35) 
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Now,  it's  a  matter  of  working  in  reverse  of  the  process  used  at  earth  departure  to  calculate 
the  delta- v.  The  burnout  velocity  at  earth  may  be  analogized  to  the  desired  reentry 
velocity  at  Mars  and  calculated  as  follows: 

v.’-v^  +  v.J  (3.36) 

The  delta-v  to  enter  a  low  Mars  orbit  is  simply  the  difference  between  this  reentry  velocity 
and  the  circular  satellite  velocity  of  the  desired  Martian  orbit: 


Avm  =  \vre  -  v(Mt/„„)i  |  (3.37) 

The  spacecraft  will  continue  to  orbit  Mars  for  as  long  as  the  astronauts  remain  on 
the  surface.  All  the  while,  both  Mars  and  earth  continue  to  revolve  about  the  sun,  so 
when  the  departure  time  arrives,  new  calculations  of  position  and  velocity  vectors  are 


necessary. 


Delta-v  at  point  B  (#2): 

The  first  and  most  important  step  in  determining  the  magnitude  of  this  third  bum  is 
to  define  the  return  orbit.  Once  again,  the  orbit  will  be  derived  from  the  time  of  flight 
constraint  imposed  on  the  return.  The  problem  consists  of  two  position  vectors  to  be 
traversed  in  a  given  time.  The  first  position  vector  is  the  location  of  Mars  at  the  time  of 
departure,  while  the  second  position  vector  is  the  location  of  the  earth  at  the  time  when 
the  spacecraft  reaches  an  earth  radius  from  the  sun.  When  considered  in  this  way,  the 
problem  may  be  identified  properly  as  a  Gaussian  problem.9  The  two  position  vectors 
accompanied  by  a  time  of  flight  exactly  define  two  closed  transfer  orbits,  one  with  a  true 
anomaly  less  than  1 80  degrees  and  one  with  an  anomaly  greater  than  1 80. 
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Subroutine  Gauss  tackles  the  problem  at  hand,  accepting  the  two  position  vectors 
and  time  of  flight,  performing  an  iterative  process,  and  returning  the  accompanying 
velocity  vectors  at  points  B  and  A.  The  vector  at  point  B  may  be  thought  of  as  v^^, 
while  the  vector  at  A  is  the  velocity  the  spacecraft  will  have  upon  earth  arrival.  Knowing 
the  heliocentric  velocity  vector  at  Mars  necessary  for  the  return  orbit  enables  you  to  treat 
this  third  burn  just  as  you  would  treat  the  original  burn  at  earth. 

The  total  deha-v  to  be  supplied  by  the  propulsion  system  is  amply  the  sum  of  the 
changes  necessary  at  these  three  points  and  represents  the  cost  of  the  time  of  flight  for  the 
mission. 


33  -  SPECIFIC  MASS  FORMULA: 


Figure  3.2  —  Block  diagram  of  NEPTUNE  loop 
The  purpose  of  the  above  sections  was  to  create  a  picture  of  what  NEPTUNE 
must  look  like,  based  on  the  mission  it  is  to  perform.  This  section  will  help  to  define  what 
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NEPTUNE  will  look  like,  based  upon  its  components.  It  will  attempt  to  define  the  whole 
as  a  sum  of  its  parts.  The  NEPTUNE  flow  is  represented  by  the  block  diagram  that 
follows. 

The  total  mass  of  the  propulsion  system  is  no  more  than  the  mass  of  its  parts.  For 
our  purposes  these  parts  include: 

H*^“HoTO  +  Htadd  +  Meenv+  Mpwrproc  +  +  H»d  (3-38) 

What  we  are  therefore  dealing  with  is  the  problem  of  assigning  values  to  this  symbolic 
formula.  Obviously,  the  masses  will  be  different  and  dependent  upon  how  demanding  the 
propulsion  requirements  are.  In  order  to  standardize  the  elements,  we  work  in  specific 
masses  so  that  the  new  formula  becomes: 

=  (3  +  3  JP.  +  PcPc  +  iUV  +  W\  +  fM>r  (3.39) 

Considering  the  NEPTUNE  mass  equation,  some  variables  are  taken  from 
previous  analysis.  The  source  power  is  defined  based  on  maximum  NERVA  capabilities. 
The  power  that  will  arrive  at  the  converter  and  power  processor  is  the  source  power, 
minus  that  power  lost  to  inefficiencies  in  the  system.  This  is  considered  in  the  efficiency 
parameter,  c.  Likewise,  the  power  that  will  arrive  at  the  thruster  is  the 
converter/processor  power,  reduced  for  inefficiencies  in  those  subsystems.  Finally,  the 
power  that  reaches  the  radiator  is  that  waste  power  lost  to  inefficiency  throughout  the 
system.  It  is  the  "other"  power. 

The  NEPTUNE  mass  equation  may  be  broken  down  to  incorporate  the  efficiency 
terms.  Consider  the  converter  subsystem,  changing  heat  from  the  source  to  energy  for  an 
electric  system.  The  amount  of  mass  this  will  contribute  to  NEPTUNE  will  look  like: 

MLr-PA  +  O-OflUP. 


(3.40) 
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ec :  converter  efficiency 

A  similar  formula  can  be  written  for  the  power  processor,  which  would  transform  the 
power  to  the  exact  form  needed  for  the  load,  the  thruster.  However,  this  is  where  one  of 
the  possible  advantages  to  the  NEPTUNE  system  comes  into  play.  It  has  been  explained 
that  NEPTUNE  derives  efficiency  through  optimum  component  coupling.  The  MHD 
converter  may  also  be  considered  the  processor  because  it  produces  electric  power  in 
exactly  that  form  required  by  its  load,  the  MPD  thruster.  As  a  result,  the  mass 
contribution  of  the  processor  might  be  eliminated.  The  mass  addition  of  the  processor 
waste  power  to  the  radiator  would  then  also  be  avoided. 

The  thruster  subsystem  is  a  unique  situation.  While  the  radiator  is  responsible  for 
maintaining  steady-state  conditions  in  the  closed  cycle  NEPTUNE  generator  and  source 
subsystems,  it  does  not  have  to  include  losses  in  the  thruster  that  are  carried  away  by  the 
exhaust  flow,  or  radiated  by  the  thruster  itself 

The  reactor  shield  is  obviously  a  direct  factor  of  the  amount  of  power  supplied  by 
the  source  reactor.  For  the  purposes  of  this  study,  the  specific  mass  of  the  shield  will  be 
combined  with  the  source  specific  mass. 

A  specific  mass  term  for  the  radiator  may  be  derived  from  an  equation  for  the  area 
of  the  radiator.  The  area  required  to  radiate  any  amount  of  power  at  a  fixed  temperature 
is  given  by  the  equation10: 

(3*1) 

T^, :  temperature  at  which  power  is  radiated 
P^  :  power  to  be  radiated 
o :  Boltzmann  constant 
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Obviously,  as  more  power  must  be  radiated,  the  area  of  the  radiator  must  increase  to 
accommodate  it.  However,  the  equation  illustrates  the  tremendous  importance 
temperature  has  on  radiator  area,  and  therefore,  radiator  mass.  The  fourth  power 
dependence  upon  temperature  means  that  the  radiating  temperature  must  be  kept  to  a 
minimum  at  all  costs  to  avoid  a  mass  explosion. 

From  the  radiator  area  and  a  knowledge  of  the  general  properties  of  the  material 

from  which  the  radiator  is  to  be  made,  a  radiator  mass  can  be  related  to  the  area: 

—  \lArad  0-42) 

p  :  specific  mass  of  radiator  (kg/mz) 

Substituting  the  equation  for  radiator  area,  a  true  specific  mass  is: 

=  £  (3  43) 

In  order  to  develop  a  single  formula  for  the  overall  mass  of  the  propulsion  system 
of  NEPTUNE,  all  of  the  terms  must  be  combined.  This  equation  can  be  reduced  so  that 
all  powers  are  defined  in  terms  of  the  known  source  power: 

Mneptuni  =  [$NERVA  +  figen&gen  +  P(6j«i6(  +  (1  ~  (3-44) 

The  above  equation  is  very  useful.  The  mass  resulting  from  the  equation  for  may 

be  compared  to  the  results  of  the  rocket  equation.  By  multiplying  the  maximum  allowable 

NEPTUNE  system  specific  mass,  from  the  mission  analysis,  by  the  NERVA  source  power, 

we  define  a  maximum  allowable  mass  for  the  NEPTUNE  system.  By  summing  actual 

specific  masses  of  the  various  components,  we  can  estimate  the  actual  mass  of  the 

NEPTUNE  system.  Ideally,  we  seek  the  inequality: 

®nepti«ie  ^  P  prep  (3.45) 

The  results  of  this  inequality  will  derive  the  ideas  presented  in  the  second  part  of 
the  thesis.  If  the  combined  mass  of  the  components  of  NEPTUNE  falls  within  the  mass 
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budget,  the  problem  of  applying  NEPTUNE  to  a  manned  Mars  mission  is  solved. 

However,  if  the  system  mass  is  above  the  allowable  maximum,  improvements  to  the 
system  must  be  engineered  in  order  to  bring  the  mass  within  limits. 


CHAPTER  IV  -  ANALYSIS  OF  TRIP-TIME  AND  SPECIFIC  MASS 


4.1  -  DELTA-V  FOR  IMPULSIVE  MARS  MISSION: 

The  goal  of  program  MARS  TRIP  (Appendix  A)  was  to  model  all  possible  round 
trip  impulsive  Mars  trajectories.  In  order  to  cover  all  possible  planetary  geometries 
coming  from  different  combinations  of  outbound  transfer  time,  surface  stay  time,  and 
return  transfer  time,  the  program  performs  delta-v  calculations  at  small  increments  of 
each.  The  procedure  outlined  in  the  previous  chapter  was  adhered  to  completely.  The 
resulting  trajectory  model  is  easily  sufficient  to  generate  the  rough  figures  necessary  for 
this  project,  but  it  clearly  incorporates  assumptions  that  separate  it  from  reality.  These 
assumptions  are  summarized  below: 

1  —  Launch  time  is  flexible.  It  is  possible  to  wait  a  launch  opportunity  so  that 

Mars  will  arrive  at  point  B  (initial  Mars  encounter)  just  as  the  transfer  craft 

does. 

2  —  The  problem  is  analyzed  in  two  dimensions. 

3  —  A  three  body  (Sun,  Mars,  earth)  scenario  is  considered. 

4  —  The  mission  begins  from  a  400  km  earth  orbit.  The  craft  enters  a  400  km  orbit 

of  Mars. 

5  —  No  final  bum  is  included  to  enter  earth  orbit  on  return. 

6  —  Trajectory  calculations  are  performed  at  twenty  day  increments.  This  is 

assumed  sufficient  to  get  clear  picture  of  delta-v  trends. 

7  -  Payload  mass  is  set  at  20,000  kg. 
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MARSTRIP  was  set  up  to  calculate  the  delta-v  requirements  for  round  trip  Mars 
missions  at  twenty  day  increments  in  outbound  transfer  time,  Martian  surface  stay  time, 
and  return  transfer  time.  The  transfer  times  were  considered  for  periods  between  70  and 
250  days.  The  stay  times  lasted  from  10  to  1 10  days.  This  resulted  in  minimum  mission 
time  of  150  days  round  trip  and  a  maximum  of  610  days.  The  minimum  was  governed  by 
the  fact  that  transfers  of  shorter  times  than  these  lead  to  unreasonably  high  delta-v 
requirements  and  hyperbolic  trajectories  with  respect  to  the  sun.  The  maximum  is 
imposed  due  to  the  fact  that  this  is  a  manned  mission  and  space  travel  nearing  two  years 
would  be  excessively  demanding. 

Figure  4. 1  and  4.2  summarize  the  results  of  MARSTRIP  when  an  exhaust  velocity 
of  98. 1  km/s  (1^=  10,000)  is  used.  The  curves  are  the  result  of  a  series  of  arcs.  Because 
of  the  possibility  for  a  return  trip  in  either  direction  about  the  sun,  the  program  calculates 
the  delta-v  for  both  options.  Therefore,  as  the  geometry  of  the  sun-earth-Mars  system 
changes,  the  direction  for  the  "short"  return  transfer  about  the  sun  reverses.  This  leads  to 
the  discontinuities  in  the  graphs.  The  graphs,  however,  clearly  highlight  the  maximum  and 
minimum  delta-v*s  and  make  it  possible  to  couple  a  delta-v  value  with  a  trip-time  value  for 
either  the  long  or  short  option.  For  more  specific  numbers,  refer  to  Appendix  B,  which 
contains  the  actual  values  used  to  generate  these  graphs. 

Specifically,  the  delta-v  requirement  for  the  fastest  mission  (150  days)  considered 
was  42.9  km/s.  The  longest  mission,  on  the  other  hand  required  a  deha-v  of  only  29.7 
km/s.  Although  it  would  seem  that  the  deha-v  and  time  of  flight  would  have  a  direct 
inverse  relationship,  this  was  not  the  case.  As  the  missions  became  less  demanding  in 
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terms  of  transfer  time,  a  general  trend  in  decreased  delta-v  requirement  was  apparent, 
however,  the  dependence  upon  convenient  planetary  geometry  was  just  as  large  a  factor. 


The  return  trajectory  was  solved  to  consider  the  option  of  a  return  in  either 
direction  about  the  Sun.  The  "long  option"  as  it  was  called  in  the  program  favored  longer 
trajectories  and  returned  the  lowest  delta-v’s.  However,  the  "short  option"  could  also  find 
small  delta-Vs  with  the  definite  advantage  that  missions  under  this  option  were  much 
shorter.  The  optimum  missions  identified  by  MARSTRIP  are  summarized  in  the  table 
below: 


Table  4. 1  —  Minimum  delta-v  round  trip  Mars  trajectories  of  each  option 


Transfer  Time 

Outbound 

Surface  Stay 

Return 

Mission  AV 

(days) 

(days) 

(days) 

(days) 

(km/sec) 

290 

130 

10 

22.5 

470 

230 

10 

230 

11.4 

Although  the  above  trajectories  result  in  the  lowest  possible  delta-v's,  they  are  not 
necessarily  the  optimum  trajectories  for  the  mission.  In  fact,  the  surface  stay  time  is  likely 
too  short  to  warrant  a  manned  Mars  excursion.  The  complete  results  of  program 
MARSTRIP  are  given  in  the  following  graphs.  They  plot  delta-v  as  a  function  of  total 
transfer  time.  Use  the  graphs  by  first  deciding  on  a  desired  round-trip  mission  time. 

Follow  the  horizontal  axis  on  either  the  short  or  long  option  graph  to  that  time.  Then 
proceed  vertically  until  a  curved  line  is  intersected.  For  most  mission  times,  several  curves 
will  correspond.  Each  represents  a  different  combination  of  outbound  transfer,  surface 
stay  time,  and  return  transfer  time.  The  discontinuities  (straight  lines  connecting  the 
curves)  in  the  graphs  come  at  the  point  where  the  "short"  trajectory  changes  to  a  return  in 
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the  opposite  direction  about  the  sun.  This  is  a  result  of  dynamic  planetary  geometry 
during  the  outbound  transfer  and  surface  stay.  Do  not  confuse  these  lines  with  the  delta-v 
curves  matched  to  the  mission  time.  For  a  specific  mission  scenario,  made  up  of  an 
outbound  transfer,  surface  stay,  and  return  transfer,  refer  to  Appendix  B.  This  is  the 
actual  output  data  from  program  MARSTRIP  and  is  not  confused  by  the  discontinuities. 

The  true  usefulness  of  the  graph  is  found  in  the  minimums  it  highlights.  The 
minimum  delta-v*s  may  be  matched  with  their  corresponding  mission  times.  The  desta-v 
parameter  may  then  be  used  set  up  a  specific  mass  budget  in  the  manner  considered  below. 

4.2  -  SPECIFIC  MASS  BUDGET: 

Next,  the  process  outlined  in  the  previous  chapter  was  used  to  calculate  a  budget 
for  the  specific  mass  of  the  propulsion  system.  Thus  budget  represents  the  largest  value  to 
which  the  specific  mass  of  NEPTUNE  per  unit  power  may  grow  before  the  mission 
becomes  impossible.  As  explained,  the  budget  will  vary  with  a  number  of  factors,  such  as 
total  bum  time  for  the  propulsion  system,  efficiency  of  the  generator  and  thruster,  and 
power  supplied  by  the  source.  Results  are  summarized  for  likely  cases  below: 

Table  4.2  --  NEPTUNE  specific  mass  budget 


AV 

P 

source 

Ue* 

fiMHD 

ffcum 

P NEPTUNE  1 

(km/sec) 

(MW) 

(km/sec) 

(%  Uj 

22.5 

98.1 

0.5 

0.5 

0.5 

13.7 

11.4 

1,732 

98.1 

0.5 

0.5 

0.5 

74.1 

22.5 

1,732 

98.1 

0.5 

0.5 

1 

39 

11.4 

1,732 

98.1 

0.5 

0.5 

1 

159.7 

22.5 

1,732 

98.1 

0.5 

0.5 

2 

89.5 

11.4 

1,732 

98.1 

0.5 

0.5 

2 

3309 
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The  source  power  value  used  in  the  calculations  comes  from  an  understanding  of 
the  power  producing  capacity  of  the  NERVA  engine.  The  exhaust  velocity  for  the  MPD 
thruster  is  based  on  projections  of  their  abilities  to  produce  I^'s  up  to  10,000  seconds. 
Efficiencies  are  very  rough  estimates  of  the  abilities  of  MHD  and  MPD  technologies.  In 
the  case  of  the  MHD  generator,  these  values  may  be  considered  conservative,  while  in  the 
instance  of  the  MPD  thruster,  they  lean  more  toward  the  liberal  side.  General  goals 
necessary  to  make  NEPTUNE  useful  would  be  to  lose  no  more  than  half  the  power  in 
transition.  Finally,  the  bum  time  is  a  factor  more  of  the  trajectory  class  considered  than 
the  ability  of  the  NERVA  engine.  Because  NEPTUNE  is  anchored  by  a  nuclear  source 
which  is  more  tailored  to  a  closed-cycle  (as  used  on  earth)  than  open,  it  should  have  little 
problem  applying  its  power  over  the  longer  term.  After  all  this  time  period  is  still  very 
short  when  compared  to  the  time  that  nuclear  generators  are  expected  to  work  on  earth. 
For  this  project,  the  bum  times  were  limited  by  the  fact  that  the  trajectories  followed  are 
impulsive  and  the  application  of  thrust  must  therefore  be  impulsive. 

A  full  plot  of  specific  mass  budget  variance  with  the  time  of  flight  data  from 
program  MARSTRIP  is  included  in  Figure  4.3  and  4.4  (for  the  short  and  long  return 
trajectory  option).  The  graph  is  obviously  very  rough  because  of  the  discontinuities  in  the 
output  of  MARSTRIP,  but  once  again,  the  idea  is  to  identify  the  maximum  budget  and 
these  are  clearly  identifiable  in  the  plot.  The  actual  output  data  contained  in  Appendix  B  is 
useful  to  obtain  numerical  values  of  the  specific  mass  budget. 

The  tradeoff  between  bum  time  and  specific  mass  is  an  interesting  one.  The  above 
table  illustrates  the  increase  in  budget  derived  from  a  longer  thrusting  time.  Figure  4.5 


DELTA-V  (km/s) 


TIME  OF  FLIGHT  VS.  DELTA-V  (SHORT  OPTION) 


200  300  400  500 

TOF  (doys) 

Figure  4.1  -  Variance  of  Deha-v  with  Mission  Time  of  Flight  for  Short  Option 
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TIME  OF  FLIGHT  VS.  DELTA-V  (LONG  OPTION) 


200  300  400  500 


TOF  (days) 

Figure  4.2  -  Variance  of  Delta-v  with  Mission  Time  of  Flight  for  Long  Option 


Specific  Moss  (kg/MW) 
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highlights  the  linear  relationship  between  these  two  parameters.  This  graph  would  seem  to 
imply  that  the  thrusting  time  should  be  increased  as  much  as  possible.  It  must  not  be 
forgotten,  however,  that  the  trajectories  being  analyzed  are  impulsive  and,  as  such,  the 
thrusting  time  must  look  instantaneous  when  compared  to  the  time  for  the  transfer.  This 
fact  dearly  limits  the  acceptable  amount  of  time  during  which  the  propulsion  system  may 
apply  its  thrust. 


4.3  -  SPECIFIC  MASS  VALUES: 

This  section  summarizes  research  on  existing  systems  which  would  serve  as 
components  of  NEPTUNE.  MHD  generators  and  MPD  thrusters  working  with 
continuous  levels  of  power  as  significant  as  those  intended  for  NEPTUNE  are  yet  to  be 
developed.  However,  through  the  use  of  the  specific  mass,  systems  designed  with  less 
ambitious  uses  in  mind  may  be  applied  to  multimegawatt  power  levels.  Obviously,  since 
the  existing  technology  is  not  developed  for  very  high  powers,  there  is  a  great  deal  of 
uncertainty  as  to  how  the  technology  will  transfer.  This  can  be  either  a  positive  or 
negative  aspect.  It  encompasses  phenomenon,  such  as  the  onset  condition,  which  could 
make  high  power  applications  impossible.  However,  more  than  likely,  a  system  designed 
to  deal  with  multimegawatt  power  levels  will  have  a  lower  specific  mass  than  their  smaller 
rdatives.  The  mass  of  the  system  will  undoubtedly  grow,  but  not  in  a  linear  fashion  with 
the  power.  As  a  result,  the  conclusions  from  this  section  are  relativdy  predictable. 

Masses  within  the  budget  cannot  be  expected  until  the  subsystems  are  reconfigured  to  deal 
with  very  high  power  systems,  such  as  NEPTUNE. 
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Because  the  source  is  set  up  for  very  high  powers,  it  is  useful  to  start  with  a 
consideration  of  NERVA.  NERVA,  acting  in  its  capacity  as  a  nuclear  thermal  rocket,  can 
possibly  attain  specific  impulses  as  high  as  1000  seconds,  translating  to  exhaust  velocities 
of  9810  m/s.  It  does  so  by  heating  the  coolant/propellant  to  a  maximum  temperature  of 
2500  K.  In  the  past  it  has  demonstrated  an  increase  in  efficiency  with  increasing  power 
output.  Although  NERVA  never  supplied  power  above  a  couple  hundred  megawatts, 
thrust  projections  for  the  system  were  as  high  as  330  kN.n  Combining  this  with  the 
specific  impulse  of  up  to  1000  seconds  would  lead  to  a  maximum  power  of  over  1600 
MW,  which  is  well  above  the  1000  MW  demands  estimated  for  this  study. 

A  number  of  different  figures  can  be  found  for  NERVA's  mass.  They  vary 
depending  on  the  output  power,  and  subsidiary  devices  included.  However,  the  specific 
mass  of  NERVA  is  about  220  kg/MW.  *  This  factor  again  is  based  on  NERVA  tests 
conducted  at  lower  power  levels  than  expected  for  NEPTUNE. 

The  NERVA  source,  because  it  is  developed  to  operate  at  fairly  high  power  levels, 
has  a  lower  specific  mass  than  the  other  subsystems  that  make  up  NEPTUNE  and  already 
it  is  apparent  that  NEPTUNE'S  system  mass,  when  developed  as  a  combination  of  the 
component  masses  will  be  larger  than  allowed  by  the  budget  constraints.  This  will  become 
even  more  obvious  when  the  specific  masses  for  existing  MPD  thrusters  and  MHD 
generators  are  included. 

In  the  area  of  MPD  thrusters,  the  specific  mass  values  vary  greatly,  due  to  the 
different  models  created.  The  self-field  vs.  the  applied-field  thrusters  are  two  completely 
different  types  of  device,  which  react  very  differently  to  high  currents.  They  are  estimated 
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to  be  capable  of  operating  in  the  specific  impulse  range  of  1000  to  10000  seconds,  but 
seem  to  work  at  power  levels  only  between  1 00  and  1 000  kW. 12  Fortunately,  however, 
their  efficiency  has  displayed  a  trend  of  improvement  at  the  higher  power  levels  (0. 1-10 
MW).  The  highest  demonstrated  performance  for  a  self-field  MPD  thruster  was  a  specific 
impulse  of  5000  seconds  when  operating  at  1 .5  MW.  The  thruster  displayed  an  efficiency 
of  40  %.  At  30  kW,  however,  the  thruster  worked  at  an  efficiency  of  70  %  with  an  1^  of 
7000  seconds.13 

Conceptually  speaking,  MPD  thrusters  are  thought  to  be  possible  with  specific 
masses  of  0.2  kg/kW.12  This  would  translate  to  200  kg/MW  and,  once  again,  push  the 
envelope  for  NEPTUNE's  budget,  even  under  best  conditions.  This  subsystem  of 
NEPTUNE  is  by  far  the  least  understood  and  warrants  further  study.  There  is  no  readily 
apparent  reason  why  MPD  thrusters  can  not  be  used  in  a  high  power  system,  and,  if  that  is 
the  case,  their  specific  impulse  values  make  them  too  attractive  to  overlook. 

The  MHD  generator  has  been  applied  widely  to  ground  based  power  producing 
systems,  so  a  great  deal  about  the  devices  is  known.  Unfortunately,  since  system  mass  is  a 
relatively  minor  constraint  on  the  ground,  mass  values  are  still  questionable.  The  MHD 
generator  is  composed  of  many  of  the  same  parts  as  the  MPD  thruster  and  can  be 
expected  to  adopt  similar  specific  masses.  Once  again,  specific  masses  for  these  systems 
will  undoubtedly  improve  with  power  level.  In  fact,  in  the  case  of  the  MHD  generator,  h 
becomes  a  far  more  attractive  system  at  higher  power  levels.  MHD  generators  mate 
extremely  well  with  nuclear  sources  because  of  the  high  energy  fluids  they  can  generate. 
Although  characteristic  values  for  specific  masses  of  these  generators  are  normally  around 
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4  kg/kW8,  one  estimate  puts  them  as  low  as  0.45  kg/kW.14  The  efficiency  of  the 
generators  vary  wildly.  The  working  fluid  and  the  ease  with  which  it  ionizes  are  major 
concerns.  Projections  of  efficiencies  of  80  %  with  cesium  gas  and  55  %  with  xenon  are 
notable.14 

4.4  -  SOLVING  THE  PROBLEM 

Combining  only  the  best  specific  masses  of  the  generator,  thruster,  and  NERVA 
source  results  in  a  NEPTUNE  specific  mass  of  870  kg/MW.  This  already  exceeds  the 
budget  defined  earlier.  This  is  a  problem  (although  not  unexpected)  considering  the  fact 
that  the  radiator  mass  has  not  been  included  in  the  estimate  and,  because  of  the 
inefficiencies  expected  in  the  system,  this  mass  could  be  significant.  Fortunately,  the 
values  used  in  the  analysis,  as  mentioned  previously,  are  for  systems  operating  at  much 
lower  power  levels.  Improvements  in  both  system  efficiency,  leading  to  lower  radiator 
mass,  and  specific  mass  will  undoubtedly  come  with  the  high  power  systems. 

The  question  of  NEPTUNE'S  ability  to  accommodate  a  manned  Mars  mission  has 
been  answered.  In  its  most  basic  state,  using  component  specific  mass  parameters 
intended  for  low  power  applications  (and  restricting  the  bum  time  to  2  percent  or  less  of 
the  trip  time),  NEPTUNE  can  not  support  this  mission.  The  question  now  becomes  -  can 
alterations  to  NEPTUNE  be  made  which  might  allow  it  support  the  manned  Mars 
mission?  This  question  is  the  topic  of  the  remainder  of  the  thesis. 


CHAPTER  V  -  NEPTUNE  DEVELOPMENT  AREAS 


The  conclusions  from  the  previous  chapter  must  be  that  if  NEPTUNE  is  to  be 
expected  to  accommodate  a  manned  Mars  venture,  one  of  two  things  will  have  to  occur  -- 
either  the  mission  must  become  less  demanding  (in  terms  of  delta-v  requirement)  or  the 
system  must  become  more  apt  to  undertake  a  demanding  mission. 

The  idea  of  reducing  the  delta-v  requirement  for  an  impulsive  Mars  mission,  while 
maintaining  a  fixed  time  of  flight  is  a  challenging  astrodynamics  problem.  One  popular 
technique  used  to  reduce  the  delta-v  burden  on  the  propulsion  system  is  the  planetary 
swingby.  The  gravity  of  large  bodies  can  be  used  to  create  a  velocity  change  on  the 
spacecraft.  Moreover,  if  done  cleverly,  this  velocity  change  can  be  appropriate  to  the 
transfer  for  which  the  craft  is  intended.  One  scenario  for  obtaining  a  "free"  delta-v  using  a 
Martian  swingby  is  addressed  in  Appendix  C.  Additionally,  a  computer  routine  to 
calculate  the  effects  of  the  swingby  is  in  Appendix  D.  This  is  a  formidable  problem  and 
deserves  further  consideration. 

This  section  will  consider  two  developments  that  could  aid  the  NEPTUNE 
propulsion  system,  itself,  in  better  using  its  power.  The  technologies  discussed  are 
intended  to  reduce  the  mass  of  NEPTUNE  without  detracting  from  its  capacity  to  produce 
thrust.  The  first  section  considers  an  alternative  MHD  generator,  in  which  a 
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nonequilibrium  plasma  is  used  to  carry  the  flow  energy.  The  second  section  considers  the 
development  of  a  different  nuclear  rocket  source,  employing  a  particle  bed  reactor. 

5.1  -  NONEQUILIBRIUM  MHD: 

The  dimensions  of  the  MHD  device  used  in  NEPTUNE  will  be  larger  than  any 
device  previously  studied  for  space  purposes  because  of  the  magnitude  of  the  power 
requirements  for  an  MPD  thruster  intended  for  interplanetary  space  exploration. 

However,  just  as  NEPTUNE  evolved  as  an  alternative  to  the  brute  force  method  of 
continually  increasing  the  size  of  established  chemical  systems,  the  more  powerful  MHD 
generator  will  be  of  little  use  if  it  simply  scales  with  the  smaller  ones  that  preceded  it. 
Therefore,  since  it  is  an  intention  of  this  project  to  improve  rather  than  to  scale,  we  must 
exploit  characteristics  of  an  MHD  generator  that  affect  its  power  output,  that  are  more 
subtle  than  its  size. 

This  section  will  present  the  theory  behind  a  nonequilibrium  MHD  generator  and 
explain  why  it  is  able  to  drive  a  high  current,  using  minimal  energy.  The  explanation  will 
describe  the  benefits  of  the  generator  by  addressing  a  series  of  related  questions  about  the 
physics  of  a  working  fluid  as  it  interacts,  first  with  an  arc  discharge,  then  with 
electromagnetic  forces  in  the  MHD  generator. 

The  first  issue  is  to  address  the  respective  importance  of  electrons  versus  ions  in 
the  working  plasma  in  carrying  current  through  the  generator.  Two  questions  should  be 
answered.  Why  are  the  electrons  in  a  plasma  almost  completely  responsible  for  carrying 
the  current  in  the  MHD  generator  and  what  relates  the  electric  field  created  in  the 
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generator  to  the  amount  of  current  that  the  generator  can  drive?  These  questions 
introduce  the  issue  of  conductivity  of  the  working  fluid  and  the  parameters  of  the  plasma 
that  influence  it.  Next,  once  it  is  established  that  a  high  conductivity  plasma  is  necessary 
to  produce  high  current,  and  having  shown  that  it  is  the  electrons  that  are  carrying  the 
current  in  the  plasma,  how  is  it  possible  to  increase  the  conductivity  of  those  elements  of 
the  plasma  which  are  responsible  for  carrying  the  current? 

First,  consider  a  plasma  flowing  through  an  MHD  generator.  As  the  conducting 
fluid  runs  normally  to  the  applied  magnetic  field  an  electric  field  is  created: 

E total  =  E applied  +  (5xi? )  (5.1) 

u  :  flow  velocity 

Assume  the  net  electric  field  generated  has  a  strength,  E.  Then  the  force  applied  by  the 
field  on  the  electron  will  have  a  magnitude,  eE,  and  will  be  directed  in  the  same  direction 
as  the  field.  Likewise,  an  ion  in  the  plasma,  under  the  influence  of  the  same  field  will  feel 
an  equal  magnitude  force  in  the  opposite  direction  due  to  the  field,  as  it  possesses  equal, 
but  opposite  charge.  However,  since  the  electron  mass,  me,  is  much  smaller  than  the  ion 
mass,  m, ,  the  acceleration  on  the  electron,  due  to  the  electric  field,  eE/me,  will  be  much 
larger  than  the  acceleration  on  the  ion,  eE/m,. 

This  describes  the  situation  in  the  plasma  as  it  passes  through  the  MHD  generator. 
After  a  fixed  time,  t,  the  electrons  have  reached  a  velocity,  based  on  their  higher 
acceleration,  which  is  much  larger  than  that  reached  by  the  ions.  In  comparison,  the  ions 
achieve  practically  no  directed  drift  velocity  due  to  the  electric  field. 

Now,  consider  the  most  fundamental  description  of  current  density.  It  is  simply 
the  amount  of  charge  passing  through  a  unit  area  per  unit  time  -  a  charge  flux. 
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Microscopically,  h  may  be  expressed  as  a  vector  summation  of  all  the  individual  charge 
motions: 

/  =1  n,q,v,  (5.2) 

l 

n :  number  density 
q :  charge 
v :  particle  velocity 

Using  the  plasma  assumption  (that  the  electron  and  ion  number  densities  are  equal)  and 
recalling  that  the  charges  on  both  ions  and  electrons  have  the  same  magnitude,  the 
majority  of  the  current  density  will  be  due  to  the  motion  of  the  faster  particles,  the 
electrons,  in  the  direction  of  the  electric  field.  In  fact,  because  the  electrons  are  so  much 
smaller  than  the  ions,  almost  all  of  the  current  will  be  carried  by  the  free  electrons  in  the 
plasma. 

The  above  description  answers  the  first  question.  It  explains  why  the  electrons, 
rather  than  the  ions,  carry  the  current  through  the  MHD  generator.  Now,  consider  the 
relationship  between  the  net  electric  field  and  the  current  density  generated.  The 

parameter  that  relates  them  is  the  bulk  conductivity  of  the  gas: 

j  =  c(E+uxB )  (5.3) 

o :  conductivity 

This  equation,  which  is  analogous  to  Ohm's  Law,  demonstrates  the  importance  of  high 
conductivity  in  creating  current.  The  calculation  of  the  conductivity  of  the  gas  is  a 
calculation  of  the  vector-averaged  velocities  of  the  elements  which  make  up  the  gas.  This 
is  no  trivial  problem  because  of  the  many  internal  and  external  forces  that  direct  them.  In 
all,  the  mass  and  charge  of  the  particles  considered,  their  state  of  random  thermal  motion 
(temperature),  the  frequency  and  detailed  characteristics  of  their  collisions  with  themselves 
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and  with  all  other  particles  in  a  gas,  and  the  prevailing  gasdynamic  flow  all  affect  the 
migration  velocity  of  each  charged  particle  and  the  bulk  property  called  conductivity  .15 

The  motion  of  the  fluid  through  the  generator  will  provide  an  electric  field  in  the 
generator.  The  purpose  of  the  generator  is  to  drive  current.  In  the  case  of  a 
multimegawatt  MPD  thruster  load,  this  current  is  significant.  Therefore,  it  is  essential  that 
the  plasma  attain  a  high  conductivity.  Because  the  electrons  are  influenced  much  more 
heavily  than  the  ions  by  field  forces,  the  conductivity  of  a  gas  is  primarily  due  to  the 
motions  of  the  electrons  through  it.16  It  is  specifically  these  particles  that  must  have  a  high 
conductivity. 

Consider  one  of  the  characteristics  of  the  gas,  the  temperature,  and  why  it 
contributes  to  the  gas  conductivity.  Conductivity  may  also  be  described  as  the  ratio  of  the 
current  density  to  the  rate  at  which  electrons  in  a  unit  volume  gain  momentum  by  impact 
with  positive  ions: 

°  *  (5  4) 

Consider  the  denominator  of  this  ratio  for  two  plasmas,  one  at  high  temperature  and  one 
at  low,  and  determine  the  effect  of  the  temperature  of  the  gas  on  conductivity.  In  the  low 
temperature  plasma,  the  electron  will  impact  ions  and  exchange  some  amount  of 
momentum  in  the  collision.  As  the  temperature  is  raised  to  the  second  situation,  the 
energy  is  increased  in  the  gas  uniformly,  so  the  velocity  increases  in  the  electrons  more 
than  the  ions  (since  it  is  a  squared  relationship  in  energy).  Now  when  an  electron  collides 
with  an  ion,  since  momentum  is  linear  with  velocity,  a  larger  amount  of  momentum  is 
imparted  on  the  ion  and  less  momentum  is  gained  by  the  electron  in  the  collision  (where 
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momentum  is  conserved).  As  a  result,  the  denominator  of  the  conductivity  ratio  becomes 
smaller  with  higher  temperature  and  the  conductivity  becomes  larger. 

A  relationship  for  the  conductivity  of  a  fully  ionized  gas  flowing  transverse  to  a 
strong  magnetic  field  (as  is  the  case  in  the  MHD  generator)  can  be  written17: 

°  =  r»*1oVZhA)  [®Wcm]  (5.5) 

Z :  atomic  number 
In  A :  h/p0 

h  :  Debye  shielding  distance 
p0 :  close-collision  distance 

This  relationship  reflects  the  fact  that  the  conductivity  is  directly  related  to  the 
temperature  of  gas. 

Now,  the  importance  of  the  electrons  to  carry  the  current  in  the  generator  has  been 
established.  Also,  the  importance  of  high  conductivity  and  the  corresponding  high 
temperature  that  helps  to  create  it  has  been  shown.  All  that  remains  is  to  find  a  way  to 
exploit  these  qualities  of  a  plasma  to  make  the  most  efficient  generator.  How  is  it  possible 
to  create  a  plasma  where  the  temperature  of  the  electrons  and  the  associated  conductivity 
is  high  enough  so  that  they  may  carry  the  necessary  amount  of  current,  but  keep  the  bulk 
temperature  of  the  gas  low  to  avoid  the  excess  energy  input  necessary  for  uniform  heating 
and  ensure  that  it  can  be  contained  by  the  materials  in  the  generator? 

The  problem  of  creating  a  plasma  with  a  much  higher  electron  than  ion 
temperature  should  begin  with  a  discussion  of  how  a  plasma  is  created,  in  general. 
Consider  the  working  fluid  in  the  NEPTUNE  system  as  it  exits  NERVA's  nozzle.  The 
energy  added  by  the  nuclear  source  to  the  working  fluid  is  used  to  create  a  high  velocity 
flow.  As  a  result,  by  the  time  the  fluid  reaches  the  generator  subsystem,  the  heat  energy 
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has  been  converted,  for  the  most  part,  to  fluid  velocity  and  the  gas  temperature  is  too  low 
for  ionization.  The  extent  to  which  magnetic  body  forces  can  be  exerted  on  an  ionized  gas 
depends  on  the  ability  of  that  gas  to  conduct  electric  current.  Obviously,  for  the  gas  to 
carry  current,  free  charges  must  exist  in  the  gas.  Processes  that  provide  these  free  charges 
in  the  gas  are  ionization  processes. 

Ionization  occurs  when  an  amount  of  energy  is  added  to  that  atom  that  surpasses 
its  ionization  potential.  In  a  gas,  the  energy  increment  necessary  to  ionize  a  constituent 
atom  may  be  delivered  in  a  number  of  different  events.  The  atom  may  undergo  an  inelastic 
collision  with  another  particle  (an  ion,  electron,  or  another  atom)  of  sufficiently  high 
kinetic  energy  : 

Xf+A-+A*+e~+X  (5.6) 

Since  the  collision  is  inelastic,  the  colliding  particle  does  not  leave  with  as  much  kinetic 

energy  as  it  came  with.  Instead  that  energy  goes  into  the  receiving  gas  atom  to  create  the 

desired  ionization. 

A  second  event  might  involve  the  atom  absorbing  an  electromagnetic  photon  of 
adequately  high  frequency: 

hv+A->A*  +  e-  (5.7) 

Alternatively,  an  electron  may  be  forcibly  extracted  from  the  atom  by  a  strong  electric 

field,  as  may  be  the  case  in  a  radio-frequency  propagating  electromagnetic  wave  train13: 
ECA)->A++e~  (5.8) 

The  energy  increment  required  to  ionize  varies  from  element  to  element.  The 

ionization  potential  of  NEPTUNE's  working  fluid,  hydrogen,  is  given  below.’6 
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Table  5.1  —  Ionization  energies  of  NEPTUNE  fuel  possibilities 


Gas 

Resonance 

potential 

eV 

Ionization 

potential 

eV 

Probable 

Ionization 

Process 

7.0 

15.38 

H2  —  H/ 

18.0 

Hj  —  IT+H 

26.0 

H 2  —  IT+H+KE 

46.0 

H,  —  HW+KE 

He 

24.58 

He  —  He+ 

Although  several  factors  influence  the  ionization  process  in  fluids  (e  g.  pressure),  the  most 
significant  of  these  is  the  type  of  fluid  to  be  ionized  and  its  associated  ionization  energy. 
Unfortunately,  the  ionization  potential  of  hydrogen  is  higher  than  most  elements,  thus 
making  ionization  more  difficult  to  achieve. 

Fortunately,  there  are  intermediate  stages  of  ionization.  The  degree  of  ionization, 
a,  of  a  plasma  is  defined  by: 

«  =  £  (5.9) 

a  :  degree  of  ionization 
iIq  :  number  density  of  particles 
ne :  number  density  of  electrons 

It  makes  sense  that  a  plasma  with  an  abundance  of  free  charges,  a  highly  ionized  plasma, 
would  be  able  to  carry  more  current  than  one  with  a  lesser  degree  of  ionization.  To  some 
extent  this  is  true.  At  degrees  of  ionization  less  than  10'3,  the  electron-ion  collisions  are 
negligible  and  the  conductivity  increases  as  a  function  of  energy  content,  reflected  by  the 
degree  of  ionization.  However,  at  degrees  of  ionization  above  10\  the  large  coulomb 
cross-section  for  electron-ion  collisions  becomes  dominant  and,  since  it  is  assumed  that  ne 
equals  n,,  the  conductivity  is  largely  independent  of  electron  density.  Therefore,  no 
significant  gain  in  electrical  conductivity,  due  to  extra  free  charges,  is  made  beyond  the 


degree  of  ionization  of  10'3.5 
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Additionally,  there  are  ways  to  reduce  the  energy  requirement  for  ionization.  One 
of  the  most  commonly  used  is  "seeding",  or  enhancing  the  working  fluid  with  some  easily 
ionizable  component  that  leads  to  further  ionization.  The  addition  of  a  seeding  element, 
such  as  cesium,  could  take  place  just  after  the  working  fluid  has  passed  through  the  core 
of  the  reactor.  This  would  allow  ionization  to  occur  at  temperatures  as  low  as  2000  K, 
which  is  below  the  temperature  of  the  gas  coming  out  of  the  core.18  Unfortunately, 
seeding  carries  with  it  certain  drawbacks.  It  adds  definite  complexity  to  the  system  as  it 
would  require  deseeding  to  occur  prior  to  the  inlet  of  the  reactor. 

In  equilibrium,  the  degree  of  ionization  in  a  gas  is  a  function  of  the  energy  content. 
Therefore,  the  degree  of  ionization  can  be  related  to  the  gas  temperature.  This  relation  is 
the  law  of  mass  action4: 

(5.10) 

ne,n,,n, :  number  densities  of  electrons, 
ions,  and  neutrals,  respectively 
k :  Boltzmann  constant 
g,go :  electronic  partition  functions  for  ions 
and  neutrals,  respectively 
m, :  mass  of  electron 
T :  electron  temperature 
s :  ionization  potential 

The  above  equation  is  useful  in  highlighting  the  elements  that  go  into  the  ionization 
process.  Unfortunately,  the  above  form  of  the  law  of  mass  action,  complete  with  its 
electronic  partition  functions,  makes  calculations  difficult.  Consequently,  for  the 
calculations  in  this  paper  a  form  of  the  Saha  equation  to  calculate  the  degree  of  ionization, 
is  used19: 

=;2r!i+2.51og10r-6.5  (5.11) 

x  :  degree  of  ionization 
p :  fluid  pressure 
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V, :  ionization  potential 
T :  electron  temperature 

In  NEPTUNE,  the  energy  necessary  to  ionize  the  working  fluid  is  added  in  a 
concentrated  fashion  by  an  arc  discharge.  The  arc  is  created  by  applying  a  high  enough 
current  between  two  electrodes  to  break  down  the  dielectric  and  allow  a  concentrated 
flow  of  electrons.  As  the  originally  inert  working  gas  comes  out  of  NERVA  and 
encounters  the  arc  discharge,  it  is  bombarded  by  the  flowing  electrons.  The  energy 
transfer  to  the  gas  comes  in  the  form  of  collisions  between  the  electrons  in  the  arc  and  the 
components  of  the  working  gas.  Collisions  with  the  comparatively  massive  nuclei  in  the 
gas  are  elastic,  with  electrons  simply  bouncing  off  and  imparting  little  energy.  Collisions 
with  the  smaller  electrons  result  in  significant  energy  transfers,  sufficient  to  ionize  the 
fluid.  Through  this  process,  the  energy  from  the  arc  is  consumed,  almost  exclusively,  by 
the  electrons  in  the  working  fluid  and  the  temperature  of  the  electrons  increases. 

The  overall  picture  that  emerges  in  the  flowing  plasma  is  the  flow  of  two  fluids 
entangled  within  one  another.  The  first  consists  of  the  electrons,  with  a  high  temperature 
corresponding  to  the  energy  input  from  the  arc  ionizer.  The  second  consists  of  the  heavy 
particles4 

This  description  answers  the  last  of  the  questions  behind  a  nonequilibrium  MHD 
generator.  It  explains  how  the  temperature,  and  associated  conductivity,  of  the  electrons 
in  the  plasma  may  be  increased  with  minimum  energy  input.  It  was  shown  earlier  that  the 
electrons  in  the  complete  plasma  (or  the  electron  gas  in  the  two  gas  model)  are  the  devices 
used  to  carry  current.  Moreover,  the  higher  the  energy  content  (the  higher  the 
temperature)  of  the  gas,  the  higher  the  conductivity  and  the  more  current  the  gas  can 
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cany.  Therefore,  if  the  gas  flowing  through  the  MHD  generator  maintains  the  flow 
velocity  of  the  gas  as  a  whole,  and  the  electrons  in  the  gas,  which  have  soaked  up  all  of 
the  energy  added  in  the  ionization  process,  maintain  their  higher  temperature,  the 
conductivity  of  only  those  particles  which  are  necessary  to  carry  the  current  will  have  been 
increased.  In  creating  this  scenario,  we  will  have  deposited  all  the  required  energy  for 
ionization  and  conductivity  exactly  where  it  is  needed  simply  because  that  is  the  way  the 
laws  of  nature  (conservation  of  energy  and  momentum)  say  it  must  happen.  This 
substantially  increases  the  efficiency  of  the  system.  All  the  benefits  of  a  high  temperature, 
high  conductivity  gas  are  being  achieved  without  the  penalty  of  a  high  energy  input. 

Unfortunately,  another  question  is  introduced.  In  order  to  achieve  sufficient 
power  output  in  the  generator,  it  is  important  that  the  nonequilibrium  condition  be 
sustained  throughout  the  length  of  the  generator.  How  long  will  this  nonequilibrium 
condition,  created  by  the  arc  discharge,  exist?  Will  there  not  be  a  tendency  for  the  fluid  to 
move  towards  equilibrium? 

It  is  here  that  a  complete  understanding  of  the  relationship  between  collision  rates 
and  energy  transfer  per  collision  becomes  important.  Using  the  arc,  an  inert  working  fluid 
which  may  be  described  as  the  flow  of  two  entangled  gases  at  a  common  temperature,  one 
consisting  of  slow  moving  nuclei  and  the  other  made  up  by  faster  bound  electrons,  has 
been  separated  and  moved  to  a  condition  where  the  temperature  of  the  heavy  gas  is 
essentially  the  same,  while  the  energy  of  the  electron  gas  is  much  higher.  The  gas  is  now  a 
nonequilibrium  gas.  The  electrons,  with  their  freshly  added  energy,  must  undergo 
collisions  between  themselves  to  adjust  to  their  new  temperature  in  a  Boltzmann 
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distribution.  (In  actuality,  the  heavy  gas  must  also  undergo  collisions  between  its  members 
to  achieve  a  higher  new  temperature  distribution,  as  it  will  soak  up  a  minute  amount  of 
energy  directly  from  the  field  as  well.)  However,  at  the  same  time  the  electrons  are 
colliding  with  other  electrons,  there  are  also  collisions  occurring  between  the  components 
of  the  two  gases,  electrons  and  ions,  tending  to  equalize  the  temperatures  between  the 
particles. 

In  order  to  determine  if  it  is  possible  to  sustain  a  nonequilibrium  gas  through  the 
generator,  the  competing  effects  which  work  to  preserve  or  eliminate  the  nonequilibrium 
must  be  examined. 


This  problem  may  be  addressed  by  analyzing  the  flow  of  energy  through  the 
electron  gas.  The  nonequilibrium  situation  is  defined  by  a  high  energy  electron  gas  and  a 
lower  energy  ion  gas.  As  energy  is  added  to  the  electron  gas,  the  nonequilibrium  becomes 
more  pronounced.  However,  as  energy  flows  from  the  electrons  to  the  rest  of  the  gas, 
equilibrium  is  approached. 

First,  consider  the  energy  leaving  the  electron  gas.  The  amount  of  power  flowing 
from  the  gas  is  a  function  of  both  collision  frequency  and  energy  transferred  per  collision. 
The  electrons  will  have  collisions  between  themselves  that  force  the  electron  gas  to  a 
Boltzmann  distribution  for  its  new  energy  content.  The  self-collision  time  between  like 


particles  may  be  written  as: 

_  \\MmTn 
,c  itZ*h>  A 


(512) 


A :  atomic  mass  number 


Intuitively,  the  terms  in  this  expression  make  sense.  As  temperature  increases,  the  energy 


and  velocity  of  the  particles  will  increase,  making  collision  more  likely  and  decreasing  the 
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time  between  collisions.  Also,  notice  the  number  density  of  the  particles  is  in  the 
denominator.  Likewise,  this  makes  sense  since  a  dense  concentration  of  particles  will  lead 
to  more  frequent  collisions  and  a  shorter  time  between  collisions.  It  is  interesting  to  note 
that,  even  in  an  equilibrium  gas  where  the  electrons  and  ions  share  a  common  temperature, 
the  time  between  collisions  for  the  electrons  will  be  shorter  than  the  heavy  particles 
because  of  the  mm  term  in  the  numerator.  As  a  result,  if  energy  is  input  uniformly  in  the 
gas,  the  electrons  will  reach  a  Boltzmann  energy  distribution  faster  than  the  heavy  particles 
since  the  amount  of  energy  per  collision  will  be  directly  scaled  in  collision  between  like 
particles.17 

However,  the  collisions  between  electrons  in  the  gas  ref  resent  an  energy  transfer 
within  the  electron  system.  The  intent  is  to  consider  the  flow  of  energy  out  of  the  system. 
This  will  come  as  a  result  of  collisions  between  the  high  energy  electrons  and  the  dower 
heavier  ions.  The  amount  of  energy  per  unit  time  being  transferred  in  these  collisions  is 
the  power  leaving  the  electron  gas  and  working  to  erase  the  nonequilibrium. 

Spitzer  defines  the  relationship: 


<rr  tx-t 

dt  t*q 


(5.13) 


This  equation  describes  how  fast  a  two  temperature  gas  equilibrates  and  assumes  a 

common  temperature.  The  result  is  the  change  in  temperature  per  unit  time.  This  can  be 

related  to  the  change  in  energy  per  unit  time  per  unit  volume  between  the  electron  gas  and 

the  heavy  particle  gas  by  substituting  the  known  relationship  for  heat  energy: 

Poutlvol  =  (5  14) 


The  result  of  the  above  equation  is  called  P^vol.  This  represents  the  flow  of  energy  out 
of  the  electron  gas.  As  it  turns  out,  the  energy  is  not  leaving  the  gas,  but  to  the  electrons, 


69 


it  looks  as  if  it  is.  Notice,  the  parameter,  t^.  This  is  the  characteristic  time  for 
equipartition  of  energy  between  the  gases.  It  is  defined  by  Spitzer  as: 


.  5.87 A.A  ( 7W  ,  £\ 3/2 

tq  n*-(Zj-ZJ)to  A  v^»-  A) 


(5.15) 

A,A^  :  atomic  mass  numbers  of  heavies 
electrons 

Z,Z^  :  atomic  numbers  of  heavies  and  electrons 
T,T^  :  temperatures  of  heavies  and  electrons 
n^. :  electron  number  density 

The  above  describes  the  collisional  process  by  which  energy  flows  from  the 
electron  gas.  Obviously,  if  this  were  the  only  input  to  the  system,  the  gas  would  work 
towards  equilibrium  and  the  usefulness  of  the  MHD  generator  would  diminish. 
Fortunately,  there  is  a  second  factor  to  consider. 

Although  inefficiency  is  usually  a  negative  side  effect  of  electric  power  generation, 
it  is  a  fortunate  side  effect  in  the  case  of  a  nonequilibrium  MHD  generator.  Refer  to  the 

MHD  generator  power  density  equation: 

Po/voI=juB-Ji  (2.7) 

The  second  term  in  this  equation  represents  the  losses  due  to  resistive  heating.  As  the 

current  is  generated  through  the  plasma,  some  of  the  energy  is  lost  due  to  resistance 


(collisions  between  electrons)  within  the  plasma.  As  a  result,  heat  energy  shows  up  in  the 
electrons  for  the  same  reason  it  went  to  the  electrons  in  the  arc  discharge.  Normally,  this 


would  be  considered  a  loss,  since  this  energy  will  be  unavailable  to  supply  the  load. 
However,  in  this  unique  case,  this  deposition  of  heat  energy  in  the  electrons  helps  to 
enforce  the  nonequilibrium  situation  that  exists  in  the  plasma.  The  input  power  per  unit 

volume  to  the  electron  gas  is  therefore: 

PJvol  =  Ji  (5.16) 

The  conductivity  of  a  plasma  is  a  complicated  variable.  Fortunately,  Spitzer  has 
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developed  a  relation  for  the  conductivity  of  a  plasma  flowing  transverse  to  a  strong 
magnetic  field: 


a=2^lpn 


(5.17) 


The  temperature  in  the  above  equation  will  be  the  electron  temperature  in  a 
nonequilibrium  gas  because  the  electrons  are  the  dements  which  carry  the  current. 

In  order  to  maintain  or  improve  the  nonequilibrium  state  in  the  plasma,  the  flow  of 
tnergy  into  the  electrons  must  exceed  the  flow  from  the  electrons  to  the  heavy  particles. 
Tins  sets  up  an  inequality  that  can  be  rearranged  to  help  identify  those  characteristics  that 


aid  in  the  maintenance  of  a  nonequilibrium  gas  in  the  MHD  generator: 


(lt«-)*r.-  <  (3.04xl0-*H 


J2 
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kZl-ZtA^) 


1/2 


(5  i?) 


The  right  hand  side  of  the  inequality  is  made  up  completely  of  known  terms  for  any  gas. 
The  left  hand  side  will  be  dependent  on  the  characteristics  of  the  generator. 

The  inequality  defines  those  characteristics  that  help  to  maintain  the 
nonequilibrium  gas.  A  low  electron  density  is  desirable.  This  makes  sense.  Fewer 
electrons  means  fewer  collisions,  leading  to  lower  energy  transfer  rate.  Additionally,  the 
lower  the  electron  temperature  the  better.  The  low  temperature  difference  leads  to  a 
smaller  energy  gradient  driving  the  gas  to  equilibrium  and  a  lower  amount  of  energy  being 
exchanged  per  collision.  Finally,  the  high  current  density  means  that  there  is  more  loss 
due  to  resistive  heating  and,  as  stated  earlier,  the  "loss"  is  a  positive  tool  for  preserving  the 
nonequilibrium. 

At  this  point,  ft  is  necessary  to  define  what  macroscopic  properties  are  conducive 
to  the  creation  of  a  nonequilibrium  gas  in  the  ionizer?  Since  ft  has  been  shown  that 
energy  added  to  the  gas  goes  almost  exclusively  to  the  electrons,  a  good  thing  is  high 


71 


energy  input  per  unit  time.  This  means  little  more  than  reasonably  high  electric  field 
strength.  Also,  it  was  shown  that  energy  is  exchanged  between  like  particles  in  the  gas 
just  as  energy  is  exchanged  between  electrons  and  the  heavy  particles  —  through  collisions. 
Therefore,  a  bad  property  is  one  that  increases  the  collision  frequency.  Temperature  is 
one  of  these  parameters,  but  to  keep  it  at  low  values  would  be  to  defeat  the  purpose  of  the 
system  as  a  whole.  However,  what  it  is  possible  to  keep  low  is  the  pressure  and  density  of 
the  flow.  A  more  loosely  packed  gas  means  that  the  distance  between  collisions  is  greater 
and  the  time  between  collisions  is  longer.  Therefore,  the  ratio  of  E/p  becomes  very 
important  for  the  creation  of  a  nonequilibrium  gas  as  can  be  seen  from  the  graph  taken 
from  Cobine.19 


Figure  5.2  ~  Creating  a  nonequilibrium  gas 


This  relationship  can  be  calculated  numerically  as  well,  by  the  equation20: 


Tt-—Tg  _  KTrtg  (X+-eE)2 

T+-  ~  (*3V)2 


(5.19) 

:  mean  free  path  length  of  electrons 
mg :  mass  of  heavy  particle 
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The  quantity  7rmg/24me  is  already  240  for  hydrogen  and  only  gets  higher  for  helium  and 
other  elements.  Added  to  that,  the  temperature  difference  gets  larger  with  the  electric 
field,  since  the  amount  of  energy  that  electrons  pick  up  per  path  length  is  XeeE  The 
influence  of  the  pressure  in  the  gas  comes  from  the  mean  free  path,  which  is  inversely 
proportional  to  pressure.  Therefore,  it  is  important  to  maximize  the  ratio  of  E/p. 

This  section  has  explained  why  a  nonequilibrium  MHD  generator  is  an 
improvement  upon  MHD  generators  which  do  not  employ  nonequilibrium  gases. 
Moreover,  it  has  identified  those  conditions  necessary  to  create  and  sustain  nonequilibrium 
conditions  in  a  MHD  generator.  The  feasibility  of  supplying  a  nonequilibrium  MHD 
generator  with  the  exhaust  gas  generated  by  the  NERVA  rocket  is  assessed  in  appendix  E. 


5.2  -  PARTICLE  BED  REACTOR: 

The  Air  Force  recently  has  supported  two  research  programs  to  produce  a  nuclear 
thermal  rocket  to  perform  missions  like  the  NERVA  rocket  of  the  late  sixties  and 
seventies.  The  first  of  these  is  a  rejuvenation  of  the  NERVA  project.  New  technologies, 
particularly  in  the  area  of  high  temperature  carbon  composites,  add  new  possibilities  for  an 
aging  system.  The  second  project  involves  the  development  of  a  compact  particle  bed 
reactor  for  space  applications.  This  is  called  the  Space  Nuclear  Thermal  Propulsion 
(SNTP)  program. 

The  NERVA  reactor  is  of  prismatic  form,  as  explained  earlier.  The  coolant  flows 
longitudinally  through  channels  surrounded  by  rods  of  fissioning  material.  Heat  is 
exchanged  from  the  fission  reaction  to  the  flowing  coolant. 
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A  particle  bed  reactor  derives  benefits  from  the  geometry  it  employs.  The  nuclear 
fuel  is  contained  in  small  spherical  shells,  coated  to  withstand  extremely  high 
temperatures.  The  fuel  spheres  are  dispersed  in  a  porous  bed  through  which  coolant  may 
flow. 

The  hydrogen  enters  the  bed  through  a  low  temperature  cold  frit.  It  picks  up  heat, 
the  byproduct  of  fission  occurring  in  the  fuel  pebbles  in  the  supporting  bed,  and  is  expelled 
through  a  very  high  temperature  hot  frit. 


Figure  5.3  -  Particle  Bed  Reactor3 

Because  the  spheres  expose  the  maximum  amount  of  area  to  the  coolant,  more  heat  can  be 
transferred  to  the  working  gas  than  in  a  more  conventional  prismatic  reactor. 

Additionally,  in  the  particle  bed  design,  the  moderator  is  separated  from  the  fuel  and 
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surrounds  the  bed  along  with  a  reflector.  Its  absence  from  direct  contact  with  the  fuel 
allows  more  direct  contact  of  the  hydrogen  with  the  fuel. 

Obviously,  an  important  consideration  in  the  development  of  this  kind  of  reactor  is 
the  ability  to  produce  materials  which  can  be  shaped  as  necessary,  but  at  the  same  time, 
withstand  temperatures  as  high  as  3000  K.  Research  in  carbon-carbon  materials  to  replace 
graphite  strengthen  the  concept.  Although  graphite  can  handle  the  high  temperatures 
inherent  to  nuclear  reactors,  the  carbon-carbon  is  a  woven  material  that,  if  engineered 
properly,  can  be  far  stronger. 

The  result  of  these  improvements  is  a  more  efficient  reactor.  SNTP  is  projected 
to  produce  coolant  temperatures  as  high  as  3000  K,  specific  impulses  somewhat  higher 
than  NERVA  (1000  seconds),  and  thrust  of 353  kN,  which  is  above  NERVA's  best 
projections.  For  our  purposes,  the  increase  in  efficiency  leads  directly  to  a  savings  in 
mass.  The  mass  of  a  particle  bed  reactor  which  can  produce  the  same  power  as  the 
NERVA  reactor  will  undoubtedly  be  lower.  Specific  mass  values  identified  for  the 
modernized  nuclear  rocket  sources  are  as  low  as  4.525  kg/MW.2' 

The  advantages  of  the  SNTP  rocket  are  an  improved  reactor  geometry  and 
utilization  of  advanced  materials.  The  drawback,  however,  is  that  this  is  not  existing 
technology.  Therefore,  the  parameters  presented  on  SNTP  are  only  projections.  Much  of 
the  project  is  still  classified  and  unavailable  for  presentation. 

SNTP  is  another  example  of  the  benefits  of  creative  engineering.  It  does  not 
involve  the  discovery  of  some  new  energy  source.  It  simply  puts  the  available  nuclear 
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energy  source  to  better  use  through  its  ideal  geometry.  This  is  an  underlying  tone  of  the 
entire  NEPTUNE  system. 


CHAPTER  VI  -  ANALYSIS  OF  NEPTUNE  SYSTEM 


The  principles  that  make  improvements  on  the  NEPTUNE  cycle  possible  were 
described  in  the  previous  chapter.  Because  the  technologies  involved,  such  as  the  particle 
bed  reactor  and  nonequilibrium  generator,  are  still  in  the  theoretical  stage  of  development, 
the  success  which  they  might  have  remains  uncertain.  Methods  may  be  created,  however, 
to  analyze  and  evaluate  the  theory  behind  these  technologies  and  determine  if  there  is  a 
possibility  for  them  to  work  as  intended.  Appendix  E  contains  an  example  of  a  theoretical 
analysis  of  the  viability  of  a  nonequilibrium  generator  supplied  by  standard  NERVA 
exhaust. 

In  this  chapter,  however,  the  assumption  is  made  that  these  improvements  may  be 
successfully  integrated  into  the  NEPTUNE  cycle.  The  intent  of  this  chapter  is  to  analyze 
NEPTUNE  by  focusing  on  the  gasdynamic  flow. 

6.1  -  THERMODYNAMIC  ANALYSIS  OF  THE  NEPTUNE  LOOP: 

NEPTUNE  is  about  energy  and  the  working  fluid  that  carries  it.  It  is  about  how 
energy  shows  up  in  the  nuclear  reactor  and  ends  up  in  the  electric  thruster.  The  picture  of 
a  flowing  gas  is  drawn  from  the  definition  of  the  gases'  properties  at  key  points  within  the 
loop.  But  what  are  the  key  points?  The  diagram  of  the  NEPTUNE  cycle  is  now  labeled 
to  identify  the  points  of  interest,  between  each  of  its  components.  The  specifics  of  the 
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components  themselves  are  not  of  primary  interest.  Instead,  the  intent  is  to  describe  the 
state  of  the  gas  after  being  acted  on  by  the  devices.  These  are  the  points  where  the 
complete  state  of  the  working  fluid  must  be  defined. 


Figure  6. 1  —  Points  of  thermodynamic  evaluation  of  NEPTUNE  cycle 


NEPTUNE'S  fluid  cycle  can  be  compared  to  a  Carnot  cycle.  A  Carnot  cycle  is  a 
fluid  cycle  consisting  of  two  isothermal  processes  and  two  adiabatic  processes.  In  the  PV 
diagram,  the  work  extracted  from  the  fluid  is  represented  by  the  area  enclosed  by  the 


curves. 
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Figure  6.2  —  Carnot  cycle22 


Obviously,  in  order  to  extract  a  large  amount  of  work,  the  temperature  change 
over  the  highest  and  lowest  point  in  the  cycle  should  be  large  —  all  of  the  energy  available 
in  the  fluid  should  be  used.  This  idea  is  captured  in  the  Carnot  efficiency,  which 
determines  how  much  of  the  available  energy  is  extracted  by  the  cycle22: 

=  (61) 

•  ¥ 

Modeling  NEPTUNE  by  a  thermodynamic  cycle  requires  some  variation  from  the 
Carnot  cycle.  The  addition  of  heat  energy  in  the  nuclear  source  is  practically  isobaric, 
while  the  stagnation  and  static  temperatures  change.  Consequently,  the  NEPTUNE  cycle 
diverges  from  the  Carnot  cycle  between  thermodynamic  points  one  and  two  (from  the 
labeled  NEPTUNE  diagram).  Next,  comes  an  adiabatic  expansion  (like  the  Carnot  cycle) 
through  the  nozzle  on  NERVA  to  exchange  the  heat  energy  for  kinetic  energy.  Following 
this,  the  gas  flows  through  the  MHD  generator  where  the  conditions  fit  into  none  of  the 
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simple  categories.  They  are  neither  adiabatic  nor  isothermal.  This  section  of  the  cycle  is 
difficult  to  define  and  will  be  considered  in  depth  later. 

After  the  generator,  the  fluid  must  be  brought  to  the  radiator  temperature  by  an 
adiabatic  diffuser  (points  5  to  6)  and  then  the  waste  heat  may  be  rejected  at  a  constant 
temperature  by  an  isothermal  radiator.  Ideally,  after  the  radiator  has  done  its  job  in 
rejecting  losses  in  the  system,  the  fluid  will  have  been  returned  to  its  initial  inlet  conditions 
for  NERVA  at  point  1 . 


Figure  6.3  -  PV  diagram  of  NEPTUNE  loop 
Some  characteristics  of  the  NEPTUNE  system  further  constrain  its  representation 
by  a  PV  diagram.  For  instance,  there  are  definite  advantages  to  maintaining  a  reactor  inlet 
temperature  (point  1)  as  high  as  possible.  Radiator  size,  and  therefore  radiator  mass,  is 
largely  dependent  upon  the  temperature  at  which  the  energy  is  to  be  radiated: 

A  _ 


(6.2) 
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Although  it  would  be  possible  to  achieve  a  higher  cycle  efficiency  with  a  larger 
temperature  spread,  it  is  far  more  detrimental,  due  to  the  fourth  power  dependence  of 
radiator  area,  to  lower  the  minimum  temperature  any  further  than  the  reactor  inlet 
requires. 

This  chapter  seeks  to  explore  an  ideal  NEPTUNE  cycle  in  which  the  radiator 
temperature  is  held  at  the  maximum  value  possible.  This  would  be  a  case  where  the 
radiator  temperature  is  equal  to  the  reactor  inlet  temperature,  requiring  no  adiabatic 
compression  of  the  fluid  to  return  it  to  reactor  inlet  conditions  (i.e.  pressure  and  density). 
Obviously,  even  this  would  be  unacceptable  unless  the  NERVA  engine  was  modified  for 
the  closed  cycle  to  allow  a  high  enough  inlet  temperature  to  keep  the  radiator  area  at 
acceptable  levels. 

Normally,  the  coolant  for  NERVA  comes  from  a  cryogenic  storage  tank,  where  it 
is  stored  at  a  very  low  temperature.  If  the  intent  of  the  analysis  is  to  complete  the  cycle 
with  a  isothermal  expulsion  of  waste  heat  by  the  radiator,  returning  the  fluid  to  its  initial 
conditions,  the  inlet  temperature  must  be  considerably  higher  than  normal.  Only  then  will 
it  permit  the  radiator  to  radiate  its  energy  at  an  acceptably  high  temperature  —  one  that 
will  keep  the  radiator  area  and  mass  within  reasonable  limits  (which  are  defined  by  the 
specific  mass  budget  derived  in  the  first  part  of  the  paper). 

Consider  then,  for  this  analysis,  a  modified  nuclear  rocket,  with  an  inlet 
temperature  of  1000  K  and  an  outlet  temperature  of 3000  K.  The  reduced  temperature 
spread  does  not  prohibit  the  rocket  from  supporting  the  MPD  load,  as  its  extremes  still 
surround  the  temperature  for  dissociation  of  hydrogen.  Therefore,  it  can  still  store  and 
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reclaim  the  energy  involved  in  dissociation  of  the  working  fluid  as  the  cycle  proceeds, 

without  making  a  large  spread  in  the  stagnation  temperature  of  the  fluid  necessary. 

The  reduced  change  in  stagnation  enthalpy  through  the  nuclear  reactor  creates  a 

need  for  a  somewhat  higher  mass  flow  rate  for  the  coolant,  in  order  to  take  out  the  heat 

energy.  The  original  mass  flow  rate  for  the  rocket  may  be  calculated  based  on  the  known 

thrust  force  and  exhaust  velocities  it  provides. 

Consider  the  SNTP  rocket  producing  a  thrust  force  of  353  kN.21  The  maximum 

specific  impulse  projection  for  the  system  is  1000  seconds,  translating  to  an  exhaust 

velocity  of  9810  m/s.  Therefore,  using  the  equation: 

Ft  -m  uei  (6.3) 

the  mass  flow  rate  may  be  calculated  as  36  kg/s.  However,  for  the  modified  rocket,  with 

the  reduced  stagnation  enthalpy  spread,  the  mass  flow  must  increase  to  50  kg/s. 

Additionally,  approximate  nozzle  dimensions  for  NERVA's  nozzle  include  a  throat 

diameter  of  0. 19  m  and  an  exit  diameter  of  4.3  m,  which  is  a  nozzle  expansion  ratio  of 

500: 1  23  Other  parameters  are  specified  for  the  reactor  chamber  (thermodynamic  point  2). 

The  chamber  temperature  is  3000  K  and  the  pressure  is  6.895xl06  N/m2  (1000  psia).21 

These  numbers,  coupled  with  assumptions  made  concerning  the  NEPTUNE  cycle, 

are  enough  to  completely  define  the  thermodynamic  state  that  must  exist  at  the  entrance  to 

the  reactor  (point  1).  As  stated,  the  pressure  will  remain  constant  between  points  1  and  2. 

Therefore,  p,  equals  6.895xl06  N/m2.  Additionally,  stagnation  conditions  are  assumed  at 

the  reactor  inlet,  so  u,  -  0.  Finally,  as  stated  earlier,  the  temperature  at  the  reactor  inlet 

will  be  1000  K  for  this  study.  This  presses  the  upper  boundaries  for  material  integrity  for 

the  reactor,  but  it  is  sufficiently  high  to  maintain  a  low  radiator  mass  and  sufficiently  low 
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so  that  the  hydrogen  coolant  is  completely  recombined.  The  perfect  gas  law  may  be  used 

to  set  the  inlet  density  that  corresponds  to  the  other  conditions. 

pi=Pi/?«7’i  (6  4) 

Next,  use  a  combination  of  known  data  from  NERVA  and  physical  laws  of 

gasdynamics  to  define  the  state  of  the  gas  at  thermodynamic  point  3,  prior  to  the 

generator.  The  exhaust  velocity  for  the  rocket  is  U3.  Assume  that  it  will  remain  constant 

across  the  ionizing  source,  so  it  may  also  be  considered  u*,  with  a  maximum  value  of 

9810  m/s  (corresponding  to  a  specific  impulse  of  1000  seconds). 

The  rocket  nozzle  exit  leads  directly  to  the  generator,  so  the  cross-sectional  areas 

at  points  3  and  4  are  the  nozzle  exit  area,  14.5  m2  for  NERVA.  This  area  will  be  taken  as 

constant  through  the  generator. 

Since  the  mass  flow,  the  area,  and  the  velocity  of  the  flow  at  the  nozzle  exit  are 
known,  they  may  be  used  to  determine  the  fluid  density  at  this  point,  from  the  law  of 
conservation  of  mass: 

P3  =  tr3  -  P<  (6.5) 

The  temperature  of  the  fluid  at  point  3  comes  as  a  result  of  applying  conservation 
of  energy  to  the  nuclear  rocket  (between  points  1  and  3).  The  equation  takes  the  form: 

/,01+^  =  Cp7’3  +  d  (6.6) 

Since  the  source  power  is: 

Pwurct  =  \mu\  (6.7) 

and  is  known  for  NERVA,  T3  may  be  calculated. 

Finally,  the  perfect  gas  law  makes  it  possible  to  determine  the  pressure  of  the 
working  fluid  at  point  3: 

Pi  -  piRmTi  =Pa 


(6.8) 


83 


Rm  :  gas  constant  for  H2 

The  conditions  at  points  3  and  4  will  be  the  same  (with  the  exception  that  the 
electron  temperature  of  the  gas  might  be  different  if  the  electrons  could  be  separately 
heated  in  a  nonequilibrium  fashion.  Also,  seeding  in  this  area  may,  in  reality,  alter 
conditions  between  3  and  4.) 

Now,  a  complete  set  of  thermodynamic  conditions  (density,  temperature,  pressure, 
and  velocity)  of  the  fluid  are  known  at  the  inlet  to  the  generator.  The  next  step  is  to 
define  a  system  of  equations  to  solve  for  the  thermodynamic  conditions  at  the  generator 
outlet. 

Since  the  area  in  the  generator  remains  constant,  the  continuity  equation  may  be 
simplified  to  provide  the  first  equation: 

p4i/4  =  P5«5  (6.9) 

Next,  apply  the  law  of  conservation  of  momentum  across  the  generator.  In  this 
calculation,  however,  it  is  essential  to  include  the  contribution  of  magnetic  pressure: 
wS  +  p4«m7’4  +  ^ 

=f»~i  +  PiRH2Ti  +  ^  (6.10) 

B  :  magnetic  field 
w  :  generator  plate  width 
d  :  generator  plate  spacing 

In  the  above  equation,  the  magnetic  field  for  the  generator  has  not  yet  been  set.  This  will 
come  as  the  result  of  the  power  requirements  placed  on  the  generator.  Assuming  some 
efficiency  for  the  MHD,  so  that  net  output  power  is: 


(6.11) 
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The  equation  for  power  production  in  an  MHD  generator  can  then  be  used  to  find  the 
necessary  value  of  B4: 

JB4dtu-J4£  =  Pgen  (6.12) 

J :  generator  current 

The  amount  of  current  in  the  generator  is  set  by  the  current  requirements  of  the 
load.  The  necessary  current  is  a  function  of  the  thrust  requirement,  which  combines  the 
available  source  power,  the  mission  delta-v,  the  attainable  exhaust  velocity  of  the  thruster, 
the  payload  mass,  and  the  specific  mass  of  the  propulsion  system. 

Ft  —  ^—[(eAw/““  —  1  )(Mpay  +  P  psP  source )]  (6  13) 

Now,  the  current  requirement  for  the  thruster  load  and  the  generator  may  be  set: 

Jmpd  =  =  Jwwfozw)  (6- 14> 

For  the  290  day  round-trip  mission  to  Mars,  described  in  Chapter  4,  using  the 
magnetoplasmadynamic  thruster  with  and  of  10,000  seconds,  the  current  is  calculated 
to  be  210  kA. 

The  change  in  the  magnetic  field  across  the  generator  is  due  to  the  combined  effect 
of  the  motion  of  the  individual  charge  carriers  in  the  plasma.  As  they  move,  they  generate 
magnetic  fields  of  their  own,  which  add  to  the  applied  field: 

B4  +  *r  =  B5  (6.15) 

Another  conservation  law  that  may  be  invoked  is  the  law  of  conservation  of 
energy.  This  must  include  a  term  representing  the  energy  flowing  from  the  system  to  drive 


the  load: 
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CpT 4  + 


nurrt 

m 


(6.16) 


cp :  specific  heat 

Finally,  the  perfect  gas  law  is  used  at  the  generator  outlet  to  relate  pressure  to  temperature 
and  density: 


Ps  =  PsRhiTs  (6.17) 

Now,  there  is  a  system  of  six  equations,  6. 7-6. 8  and  6. 10-6. 13,  to  solve  for  the  six 

unknown  quantities: 

ps,  Ts,P5tUs,B4,Bs 

Solving  them  results  in  a  complete  definition  of  the  thermodynamic  conditions  at  the  outlet 
to  the  MHD  generator,  point  5. 

Having  addressed  the  unusual  situation  of  the  MHD  generator,  the  next  step  is  to 
continue,  using  variation  of  the  same  conservation  equations,  and  close  the  NEPTUNE 
thermodynamic  loop.  This  leads  first  to  an  analysis  of  the  adiabatic  difluser  that  is  used  to 
compress  the  flow  to  the  desired  radiator  temperature. 

Earlier,  the  conditions  which  help  to  define  an  acceptable  radiator  temperature 
were  presented  and  a  temperature  was  set  at  1000  K.  Using  this  fact,  the  mach  number  of 
the  flow  at  point  6,  just  prior  to  the  radiator,  can  be  defined  using  the  stagnation 
temperature  of  the  fluid: 

Me  =  (6.18) 


As  no  energy  has  been  exchanged  between  points  5  and  6,  the  stagnation  temperatures  of 
these  two  points  will  be  the  same.  Therefore,  the  thermodynamic  conditions  defined  at 
point  5  can  be  used  to  determine  this  stagnation  temperature: 


#5  =  JjRfiTs 

_  “5 


(619) 

(6.20) 
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7'os  =  r5(  1  +  JpMf) 


(621) 


Similarly,  the  stagnation  pressure  at  points  5  and  6  will  be  the  same  and  can  be  calculated: 

Pm  =/>s(l  +  (6.22) 

The  determination  of  the  stagnation  pressure  makes  it  possible  to  solve  for  the  pressure  of 

the  fluid  at  point  6: 


P  6  = 


P06 

(l+^A/yn 


(6.23) 


Finally,  the  perfect  gas  law  may  be  used  to  find  the  density  of  the  fluid  at  point  6: 


P6  = 


P  6 


(6.24) 


and,  with  the  local  speed  of  sound,  the  velocity  of  the  fluid  is  found  at  point  6  to  complete 


the  thermodynamic  definition  of  the  fluid: 

#6  =  JyRmTe  (6.25) 

U(,  =  (6.26) 

The  above  analysis  has  taken  one  more  step  along  NEPTUNE'S  thermodynamic 


cycle.  Now,  a  complete  set  of  thermodynamic  conditions  in  the  working  fluid  has  been 
defined  for  points  3  through  6.  Additionally,  the  conditions  that  must  exist  at  the  rocket 


inlet  are  known.  The  theory  is  that  an  isothermal  rejection  of  heat  by  the  radiator 
(between  point  6  and  1)  can  lead  directly  to  the  inlet  conditions  required  to  complete  the 
cycle.  Now,  it  is  time  to  test  that  theory. 


The  energy  equation  for  the  radiator  region  is  replaced  by  the  assumption  that 
T^T,.  However,  an  equation  for  the  conservation  of  momentum  is  available  and  takes  the 
form: 

Ui-Ui  =  RmT6  tag  (6.27) 

For  complete  stagnation  of  the  flow  at  point  1,  u,=0,  as  stated  earlier.  (Note  that  the  heat 
rejected  during  this  isothermal  process  is  ^  =  Rm  7«ln  g  .)  Again,  the  perfect  gas  law 
may  be  used  to  provide  the  last  unknown  state  variable: 
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Pi^PiRmh  (6.28) 

All  that  remains  is  to  interpret  the  results.  What  will  emerge  from  the  analysis  of 

the  cycle  as  it  proceeds  forward  (across  the  radiator)  is  a  definition  of  the  state  of  the  fluid 
at  the  inlet  to  the  reactor.  However,  from  a  knowledge  of  the  NERVA  outlet  and 
chamber  conditions  and  assumptions  of  the  flow  across  the  reactor,  a  complete  set  of 
thermodynamic  conditions  that  must  exist  in  the  fluid  at  point  1  is  already  known. 
Obviously,  if  the  cycle  is  to  be  closed,  these  conditions  must  agree. 

But  what  if  the  conditions  don't  match?  Then  it  is  necessary  to  consider  why  they 
have  not  matched.  What  characteristics  of  the  cycle  can  contribute  to  rectifying  the 
discontinuity? 

If  the  fluid  density  and  pressure  can  not  be  returned  to  their  initial  value,  when  the 
temperature  is  set,  then  there  must  not  be  enough  energy  remaining  to  complete  the  cycle. 
The  parameter  that  dictates  the  amount  of  energy  left  in  the  fluid  is  the  efficiency  of  the 
generator.  It  is  likely  that  the  system  will  be  too  efficient  to  work  in  a  closed  cycle.  In 
that  case,  the  efficiency  should  be  reduced  to  accommodate  the  closed  cycle. 

The  above  process  was  applied  to  a  system  with  a  generator  efficiency  of  0.5.  In 
this  case,  the  fluid  indeed  fell  short  of  returning  to  its  inlet  conditions.  To  correct  the 
problem,  different  generator  efficiencies  were  evaluated  until  a  satisfactory  result  was 
reached  at  an  efficiency  of  0.38.  This  efficiency  value  allowed  enough  energy  to  be 
retained  in  the  flow  so  that  the  calculated  inlet  conditions  to  the  system  were  closely 
matched  to  the  required  inlet  values  and  the  cycle  was  successfully  closed. 
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Table  6.1  --  Thermodynamic  data  of  NEPTUNE  loop 


Test  Point 

Mass  Flow 
(k«/s) 

Velocity 

(m/s) 

Density 

(kft/m3xl0^) 

Temperature 

(K)_  _ 

Pressure 

(N/m2) 

1 

50 

1.66 

1,000 

6. 895x10s 

2 

50 

— 

0.28 

3,000 

6.895xl06 

3 

50 

9,810 

3.51X10-4 

1,000 

1,459.7 

4 

50 

9,810 

3.51X10-4 

1,000 

(T^=  10,000) 

1,459.7 

5 

50 

5,320 

6.48x10“* 

2,416 

6,504 

6 

50 

8,320 

7.30x1  O'5 

1,000 

303.27 

CHAPTER  Vll  -  CONCLUSIONS 


This  thesis  set  out  to  analyze  the  ability  of  the  NEPTUNE  propulsion  concept  to 
support  a  manned  Mars  venture.  A  Mars  mission  is  a  formidable  task.  NEPTUNE  is  an 
electric  propulsion  concept,  acting  in  a  thrust  region  normally  reserved  for  chemical,  or 
possibly  nuclear  thermal  systems.  This  is  a  thrust  region  that  the  multimegawatt 
NEPTUNE  system  can  potentially  reach.  However,  the  analysis  is  clouded  by  the  fact  that 
the  multimegawatt  elements  that  make  up  NEPTUNE  do  not  exist.  Instead,  a  knowledge 
about  these  elements  must  be  implied  from  a  knowledge  of  their  lower  power  relatives. 

The  results  of  the  first  part  of  this  thesis  reflect  the  fact  that  NEPTUNE  can  not 
support  the  manned  Mars  mission,  particularly  if  we  simply  accept  component  mass 
estimates  based  on  direct  upscaling  of  the  masses  of  similar  low  power  components. 

(Bear  in  mind  that  this  conclusion  also  incorporates  thrusting  times  of  no  more  than  2%  of 
the  mission  time,  as  a  limitation.  There  is  evidence  that  suggests  that  thrusting  times 
lasting  as  long  as  5-10  percent  of  the  mission  time  may  still  be  considered  by  impulsive 
analysis.52)  Therefore,  in  order  to  rectify  the  situation,  one  of  two  things  must  be  done. 
The  manned  Mars  mission  must  become  less  demanding,  or  proof  must  be  provided  that 
multimegawatt  extrapolation  of  the  elements  which  make  up  NEPTUNE  involves  more 
than  the  direct  application  of  specific  mass  values  for  their  low  power  relatives. 
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Appendix  C  addresses  the  possibility  of  decreasing  the  delta-v  demand  of  a 
manned  Mars  venture,  perhaps  making  it  more  reachable  by  NEPTUNE.  The  second  part 
of  the  thesis,  however,  is  dedicated  to  increasing  the  effectiveness  of  NEPTUNE.  Two 
primary  concepts  were  discussed  —  the  nonequilibrium  MHD  generator  to  produce 
electricity  for  the  thruster,  and  a  particle  bed  nuclear  thermal  rocket  to  produce  high 
velocity  plasma  flow  for  the  generator.  The  reasons  why  each  device  should  improve  the 
NEPTUNE  system  were  presented. 

The  theory  behind  the  elements  which  would  make  a  Martian  voyage  possible 
using  NEPTUNE  was  summarized  by  a  thermodynamic  analysis  of  NEPTUNE  which 
incorporated  the  improvements.  This  analysis  was  intended  to  demonstrate  that  if  the 
theories  explained  could  be  incorporated  into  NEPTUNE,  the  system  could  act  on  a 
closed  cycle  to  meet  the  power  needs  of  its  Mars  mission. 

For  the  purpose  of  a  thermodynamic  analysis.  Chapter  6  assumes  a  successful 
nonequilibrium  MHD  generator.  Appendix  E  takes  a  closer  look  at  the  actual  feasibility  of 
a  nonequilibrium  MHD  generator  for  NEPTUNE  based  on  elevated  electron  temperature 
and  develops  a  method  to  analyze  the  fluid's  ability  to  sustain  a  nonequilibrium  state 
(T>>T,). 

It  is  nearly  impossible  to  apply  quantitative  values  defining  the  effectiveness  of  the 
advanced  technologies  without  actual  experimental  analysis.  The  intent,  however,  was  to 
describe  the  physics  behind  the  developments  so  that  a  legitimate  claim  may  be  made  —  the 
technologies  must  work  because  the  laws  of  physics  say  they  must  work.  Obviously, 
solutions  don't  come  that  easily.  The  basic  laws  of  nature  are  clouded  by  competing 
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phenomena  along  the  way.  The  experimentalist  must  maintain  a  belief  in  physics  and 
know  that  the  device  can  succeed,  so  that  he  or  she  can  develop  methods  in  which  the 
device  does  succeed. 

NEPTUNE  capitalizes  on  several  basic  ideas  which  immediately  make  it  a 
propulsion  system  worth  studying.  It  is  clearly  a  system  which  makes  an  effort  to  match 
components  that  belong  together,  because  of  their  separate  characteristics.  It  is  a  system 
that  uses  the  best  known  energy  source,  in  terms  of  power  density.  Finally,  it  is  a  system 
that  capitalizes  on  natural  electromagnetic  field  forces  to  employ  its  power.  Together, 
these  components  make  up  an  unproven  and  futuristic  system,  but  the  ideas  upon  which  it 
is  based  clearly  define  it  as  a  system  representing  the  direction  of  space  propulsion  on  the 
horizon. 


APPENDIX  A 


Computer  program  MARSTRIP  is  to  calculate  the  delta-v  requirement  for  a 
round-trip  Mars  trajectory,  based  on  an  input  time  of  flight  constraint.  It  cycles  through 
all  of  the  time  of  flight  and  planetary  geometry  scenarios  reasonable  for  a  manned  mission 
The  program  is  written  in  Fortran  and  frequent,  descriptive  comments  are 
contained  to  outline  the  flow.  The  program  incorporates  some  subroutines  from  the 
United  States  Air  Force  Academy  Department  of  Astronautical  Engineering  Computer 
Library.24 
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APPENDIX  B 


This  section  contains  the  two  output  files  generated  by  program  MARSTRIP, 
ct  mtained  in  Appendix  A.  The  scenario  is  based  on  an  input  exhaust  velocity  of  98. 1 
km/s.  The  output  files  represent  the  same  outbound  Martain  trajectory  and  differ  in  the 
direction  of  their  return  trajectory.  The  data  in  these  files  was  used  to  generate  the  graphs 
contained  in  Chapter  IV  of  the  thesis. 
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APPENDIX  C  -  FEASIBILITY  OF  SCHWING  TRAJECTORY  FOR  NEPTUNE 


There  are  two  ways  that  a  spacecraft  is  manipulated  as  it  follows  a  trajectory.  The 
most  obvious  is  through  the  use  of  its  on  board  propulsion  system.  Obviously,  this  has  its 
costs.  For  the  energy  to  be  used,  it  must  first  be  carried  in  some  state  and  processed  into 
another.  The  second  force  is  far  simpler.  It  results  from  the  body  forces  of  the  massive 
objects  in  range  which  assert  their  influence.  Swingbys  capitalize  upon  this  attraction  to 
obtain  "free"  delta-v's  or  delta-v's  of  convenience.  If  the  spacecraft  passes  close  to  the 
surface  of  a  much  larger  body,  the  body  will  have  the  effect  of  perturbing  the  trajectory  of 
the  vehicle  even  if  it  continues  to  revolve  on  the  larger  scale  around  another  larger  body. 
In  other  words,  this  body  can  exert  influence  without  actually  capturing  the  craft. 


107 


108 


Just  by  looking  at  the  geometry  of  the  problem,  it  can  be  seen  that  a  craft  passing 
in  front  of  the  direction  of  motion  of  the  larger  body  will  experience  a  delta-v  that  serves 
to  slow  its  global  velocity  with  respect  to  the  larger  central  body.  Conversely,  if  it  passes 
behind  the  path  of  the  perturbing  body,  its  velocity  with  respect  the  central  body  will 
increase.  (This  may  be  seen  by  comparing  the  magnitude  of  the  vector  v^,^,  in  the  two 
cases  shown  in  the  diagram.) 

Schwing  was  developed  to  reduce  the  delta-v  requirement  for  a  round-trip  mission 
to  Mars.  It  benefits  from  the  fact  that  the  transfer  craft  never  has  to  enter  a  Martian  orbit. 
After  following  a  fairly  high  energy  outbound  transfer,  the  spacecraft  drops  a  manned 
landing  payload  off  at  an  initial  Mars  encounter.  The  proximity  of  this  encounter  will 
vary,  but  will  be  well  within  the  Martian  sphere  of  influence.  As  shown  above,  the 
trajectory  of  the  spacecraft  is  perturbed  by  Mars.  The  key  is  to  perturb  the  trajectory 
appropriately  to  align  the  new  trajectory  with  a  Martian  reencounter  at  about  the  same 
time  the  landing  party  is  ready  for  departure. 

However,  the  savings  are  not  complete.  At  the  second  Martian  encounter,  the 
transfer  vehicle  is  conveniently  aligned  to  use  Mars  again.  The  purpose  this  time  is  to 
align  the  return  trajectory  so  that  a  smaller  delta-v  is  necessary  to  meet  earth  after  the 
return  time  passes. 

The  mathematics  involved  in  dealing  with  planetary  swingbys  are  greatly  simplified 
if  the  problem  is  broken  down  into  the  appropriate  reference  frames.  As  the  spacecraft 
approaches  Mars,  it  feels  no  attraction  from  the  planet  and  follows  a  completely 
heliocentric  trajectory.  Then  at  the  point  where  the  spacecraft  is  first  able  to  feel  the 
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gravitational  attraction  of  Mars,  call  it  r^,  we  define  the  beginning  of  a  hyperbolic 
trajectory  about  Mars.  The  velocity  at  this  point  in  the  heliocentric  trajectory  can  be 
considered  (v^,^^,  or  the  velocity  of  the  satellite  with  respect  to  the  sun.  In  order  to 
find  this  (v^^,  simply  use  the  relation: 


(^wtmxn)-oo  (^ui/iud)-oo  ’  Si 


^tnan/sun 


(Cl) 


(vmtolnl  00 :  velocity  of  craft  with  respect  to 
Mars  at  initial  point  of  influence 

The  next  step  is  to  find  the  eccentricity  of  this  hyperbolic  swingby  orbit  about 
Mars.  In  order  to  do  this,  the  minimum  swingby  altitude  must  be  known.  Obviously,  the 
closer  the  craft  gets  to  the  surface,  the  stronger  the  force  of  the  gravitational  attraction 
which  acts  to  deflect  the  trajectory  and  the  larger  the  deflection.  With  this  parameter 
defined,  the  eccentricity  is: 

^=i+-vzi-r  cc2> 

:  radius  of  periapsis 

This  eccentricity  makes  it  possible  to  solve  for  the  deflection  angle  in  the  hyperbolic  orbit 
caused  by  the  gravitational  acceleration  of  the  planet.  This  formula  for  this  calculation  is: 
sin(f)  =  7  (C3) 

Knowing  the  deflection  angle,  it  is  possible  to  determine  the  velocity  vector  of  the 
spacecraft  with  respect  to  Mars  as  it  leaves  the  Martian  sphere  of  influence.  The 
magnitude  of  this  vector  will  remain  as  it  was  at  the  point  of  entrance  into  the  sphere  of 
influence,  however  the  direction  of  the  vector  will  be  a  simple  rotation  through  the  angle  5 
from  the  direction  of  entrance.  As  it  turns  out,  the  deflection  angle  will  be  the  same 


whether  the  craft  passes  in  front  of  or  behind  the  path  of  the  planet.  The  difference  will  be 
in  the  direction  of  the  rotation  of  the  velocity  vector. 

Once  the  vector  velocity  of  the  craft  with  respect  to  Mars  is  known,  the 
heliocentric  velocity  of  the  vehicle  at  that  point  where  it  departs  the  Martian  sphere  of 
influence  is  found  with  simple  vector  addition: 

(Xm/sun^-oe  0£ut/mm)-oo  ^nn/ui 

The  key  to  using  Schwing  properly  is  to  determine  what  magnitude  and  direction 
of  velocity  is  desired  in  a  heliocentric  sense  upon  leaving  the  planetary  sphere  of  influence. 
In  the  case  of  the  first  passage,  the  desired  velocity  vector  is  one  that  places  the  satellite 
close  to  the  position  of  Mars  after  the  surface  stay  time  has  passed  in  a  time  of  flight 
relatively  close  to  that  surface  stay  time.  The  second  swingby  is  intended  to  leave  the  craft 
with  a  velocity  vector  enabling  it  to  conveniently  intercept  earth's  orbital  distance  from  the 
sun  at  that  time  when  the  earth  will  be  there. 

Ideally,  Schwing  would  make  it  necessary  to  bum  only  once  for  the  entire  round 
trip  mission  to  Mars.  After  the  bum  to  initiate  the  outbound  trajectory,  the  craft  would 
proceed  to  Mars  where  it  would  swing  into  a  new  trajectory  reintercepting  Mars  after  the 
surface  stay  is  complete.  Then,  it  would  swing  into  a  final  return  orbit  to  rendezvous 
with  the  earth.  Obviously,  the  gravitational  perturbing  effect  of  Mars  is  insufficient  to  act 
alone  in  steering  and  accelerating  the  craft  with  respect  to  the  sun,  but  it  is  extremely 
useful  in  reducing  the  amount  of  input  the  thrusters  must  have  to  complete  the  mission. 


APPENDIX  D 


Computer  program  SCHWING  is  to  calculate  the  delta-v  requirement  of  a 
round-trip  Mars  Schwing  trajectory.  It  is  set  up  to  define  the  delta-v  for  a  single  time  of 
flight  input,  rather  than  cycle  through  all  possible  time  of  flight  and  planetary  geometry 
scenarios,  as  in  program  MARSTRIP. 

This  program  is  written  in  Fortran  and  contains  frequent,  descriptive  comments  to 
outline  the  flow.  It  also  contains  some  subroutines  taken  from  the  USAFA  Department  of 
Astronautical  Engineering  Computer  Library.24 
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APPENDIX  E  -  FEASIBILITY  OF  A  NONEQUILIBRIUM  MHD  FOR  NEPTUNE 
The  ability  to  create  and  sustain  a  nonequilibrium  state  (with  electron  temperature 
»  heavy  particle  temperature)  throughout  the  MHD  generator  is  highly  dependent  upon 
the  density  of  the  gas  at  the  entrance  to  the  generator.  This  is  due  to  the  fact  that  a  more 
dense  gas  will  have  more  collisions  and  tend  to  equilibrate  more  quickly.  This  analysis  will 
determine  if  the  density  of  the  working  fluid  leaving  the  NERVA  nozzle  is  low  enough  to 
sustain  a  nonequilibrium  state  through  the  generator. 

The  characteristics  of  NERVA  described  in  Chapter  6  explained  that  for  a  nozzle 
exit  area  of  14.5  m2,  the  density  of  the  working  fluid  leaving  the  system  would  be 
3. 5  lxl O'4  kg/m3.  With  hydrogen  gas  as  the  working  fluid,  this  translates  to  a  number 
density  of  1 ,058xl017  cm3. 

Return  now  to  the  inequality  derived  in  Chapter  5.  It  explained  the  conditions  that 

would  have  to  exist  in  the  gas  in  order  to  sustain  a  nonequilibrium  state  in  the  generator: 
(w«-)27--  <  (S.OtxlO-4^ 

J2  kZiZ(A-)'*  {*l) 

Recall  that  the  right  side  is  completely  defined  by  the  working  fluid.  For  hydrogen  gas, 
the  inequality  may  be  rewritten: 

<>  7.33  x  1013  (E2) 
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Now,  apply  the  same  constraints  used  to  define  the  thermodynamic  cycle  in 
Chapter  6.  Additionally,  substitute  ne_  =  an,  and  solve  for  the  required  number  density 
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that  must  exist  in  the  gas  for  a  nonequilibrium  condition  to  be  sustained. 


J2(7.33x1013) 

rr  £ - ; — — 

(hv)2Tt-a? 


(E3) 


J  :  MHD  current 

The  current  required  to  support  the  mission  analyzed  thermodynamically  in  Chapter  6  is 


210  kA. 


For  reasons  explained  in  Chapter  5,  the  desired  degree  of  ionization  in  the  gas  is  at 
least  10'3.  The  Saha  equation  is  useful  in  defining  the  electron  temperature  necessary  to 
achieve  that  degree  of  ionization  in  a  hydrogen  working  fluid.  It  reveals  that  electron 
temperatures  in  excess  of  10,000  K  will  be  required  to  sustain  the  desired  degree  of 
ionization  at  this  density. 

The  only  parameter  left  undefined  in  the  inequality  is  the  channel  length  of  the 
generator  and  the  width  of  the  generator  plates.  First,  consider  a  case  where  the  width  of 
the  plates  and  the  plate  spacing  is  equal.  Since  the  cross-sectional  area  of  the  nozzle  inlet 
is  known  from  NERVA  (4.3  m),  the  width  is  defined  at  2.07  m.  Additionally,  consider  a 
channel  length  of  1  m.  Inputting  these  values  into  the  inequality  sets  the  maximum  number 
density  for  the  gas  at  8.29xl012  cm'3. 

Comparing  this  figure  to  the  actual  number  density  at  NERVA  exit,  it  is  apparent 
that  the  gas  is  considerably  too  dense  (four  orders  of  magnitude)  to  sustain  a 
nonequilibrium  state  in  the  generator.  However,  looking  again  at  the  inequality,  the 
geometrical  dimensions  of  the  generator  itself  can  contribute  to  its  ability  to  sustain 
nonequilibrium. 
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It  is  obviously  advantageous  to  create  a  generator  with  short  channel  length  and 
small  plate  width.  In  other  words,  the  generator  plate  area  should  remain  small,  leading  to 
a  large  current  density  and  more  resistive  heating.  In  order  to  create  a  small  plate  width, 
the  plate  spacing  must  grow  to  retain  the  same  cross-sectional  area. 

Consider  a  second  case  where  the  channel  length  is  reduced  to  0. 1  m.  Also,  let  the 
plate  width  be  a  quarter  of  the  plate  spacing,  w  =  1 .04  m.  Now,  the  maximum  number 
density  increases  to  1 ,73xl014  cm'3 .  Still,  however,  there  is  more  than  a  two  order  of 
magnitude  differential  between  the  required  and  actual  densities. 

In  order  to  illustrate  the  significance  of  this  differential,  an  exit  area  requirement  on 
the  nozzle  of  NERVA  which  would  result  in  the  desired  density,  may  be  calculated. 

Since  the  working  fluid  may  be  considered  completely  recombined  into  at  the  original 
nozzle  exit,  the  equations  for  an  adiabatic  expansion  through  a  nozzle  apply  beyond  this 
point.  Using  the  law  of  conservation  of  mass  and  holding  the  velocity  at  the  nozzle  exit 
constant  at  9810  m/s: 

A  required  =  mmnnqiniu  (E4) 

The  area  of  a  nozzle  which  produces  flow  rarefied  enough  to  sustain  nonequilibrium  in  the 
MHD  generator  is  6390  m2.  This  is  an  outlet  radius  of  45  meters. 

This  is  obviously  an  unfortunate  result,  as  it  would  represent  tremendous  savings 
to  be  able  to  use  a  nonequilibrium  MHD  generator  in  NEPTUNE.  However,  an 
alternative  method  of  achieving  the  necessary  conductivity  in  the  working  fluid  would  be 
to  seed  the  fluid  with  liquid  metal  droplets  (e.g.  lithium).  This  can  increase  the 
conductivity  of  the  working  fluid  without  making  extreme  temperatures  necessary. 
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Additionally,  different  generator  geometry  could  be  used.  Creating  a  series  of 
smaller  generators,  combining  to  support  the  load,  is  a  possibility.  The  small  generators 
would  have  high  current  densities  and  more  power,  via  resistive  heating,  would  be  input  to 
the  working  fluid. 


Figure  El  —  Series  of  MHD  Generators 
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